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ABSTRACT 

An  abstract  framework  for  the  analysis  and  approximation  of  a  class  of  nonlinear  optimal  control 
and  optimization  problems  is  constructed.  Nonlinearities  occur  in  both  the  objective  functional  and 
in  the  constraints.  The  framework  includes  an  abstract  non!  .ear  optimization  problem  posed  on 
infinite  dimensional  spaces,  an  approximate  problem  posed  on  finite  dimensional  spaces,  together 
with  a  number  of  hypotheses  concerning  the  two  problems.  The  framework  is  used  to  show  that 
optimal  solutions  exist,  to  show  that  Lagrange  multipliers  may  be  used  to  enforce  the  constraints, 
to  derive  an  optimality  system  from  which  optimal  states  and  controls  may  be  deduced,  and  to 
derive  existence  results  and  error  estimates  for  solutions  of  the  approximate  problem.  The  abstract 
framework  and  the  results  derived  from  that  framework  are  then  applied  to  three  concrete  control 
or  optimization  problems  and  their  approximation  by  finite  element  methods.  The  first  involves 
the  von  Karman  plate  equations  of  nonlinear  elasticity,  the  second  the  Ginzburg-Landau  equations 
of  superconductivity,  and  the  third  the  Navier-Stokes  equations  for  incompressible,  viscous  flows. 
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1.  Introduction 


The  need  to  solve  optimization  and  control  problems  arises  in  many  settings. 
Although  in  some  cases  these  problems  may  be  easily  solved,  either  analytically  or 
computationally,  in  many  other  cases  substantial  difficulties  are  encountered.  For 
example,  candidate  optimal  states  and  controls  may  belong  to  infinite  dimensional 
function  spaces  and  one  may  have  to  minimize  a  nonlinear  functional  of  the  state  and 
control  variables  subject  to  nonlinear  constraints  that  take  the  form  of  a  system  of 
partial  differential  equations  whose  solutions  are  in  general  not  unique.  In  this  paper, 
our  goal  is  to  construct  and  analyze  a  framework  for  the  approximate  solution  of 
many  such  problems.  The  setting  for  our  framework  is  a  class  of  nonlinear  control  or 
optimization  problems  which  is  general  enough  to  be  of  use  in  numerous  applications. 
The  major  steps  in  the  development  and  analysis  of  our  framework  are  as  follows: 

•  define  an  abstract  class  of  nonlinear  control  or  optimization 
problems; 

•  show  that,  under  certain  assumptions,  optimal  solutions  exist; 

•  show  that,  under  certain  additional  assumptions,  Lagrange  mul¬ 
tipliers  exist  that  may  be  used  to  enforce  the  constraints; 

•  use  the  Lagrange  multiplier  technique  to  derive  an  optimality 
system  from  which  optimal  states  and  controls  may  be  deduced; 

•  define  algorithms  for  the  approximation,  in  finite  dimensional 
spaces,  of  optimal  states  and  controls;  and 

•  derive  estimates  for  the  error  in  the  approximations  to  the  op¬ 
timal  states  and  controls. 

Two  of  the  key  ingredients  used  to  carry  out  the  above  plan  are  a  theory  given  in 
[21]  for  showing  the  existence  of  Lagrange  multipliers  and  a  theory  first  developed  in 
[6]  for  the  approximation  of  a  class  of  nonlinear  problems.  In  both  of  these  theories, 
certain  properties  of  compact  operators  on  Banach  spaces  play  a  central  role.  We 
point  out  that  the  nonuniqueness  of  solutions  of  the  nonlinear  constraint  equations 
deems  it  appropriate  to  employ  Lagrange  multiplier  principles. 

After  having  developed  and  analyzed  the  abstract  framework,  we  will  apply  it 
to  some  specific,  concrete  problems.  In  each  case,  we  use  the  abstract  framework  to 
analyze  the  concrete  problems  by  merely  showing  that  the  latter  fit  into  the  former. 
The  particular  applications  we  consider  are: 

•  control  problems  in  structural  mechanics  having  geometric  non- 
linearities  that  are  governed  by  the  von  Karman  equations; 

•  control  problems  in  superconductivity  that  are  governed  by  the 
Ginzburg- Landau  equations;  and 

•  control  problems  for  incompressible,  viscous  flows  that  are  gov¬ 
erned  by  the  Navier-Stokes  equations. 

In  considering  these  applications,  we  will  purposely  choose  different  types  of  controls 
in  order  to  illustrate  how  these  can  be  accounted  for  within  the  abstract  framework. 
In  all  three  cases,  approximation  will  be  effected  through  the  use  of  finite  element 
methods. 
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2.  The  abstract  problem  and  its  analysis 


In  this  section  we  define  and  analyze  an  abstract  class  of  constrained  nonlinear 
control  problems;  an  outline  of  the  definitions  and  results  of  this  section  is  as  follows. 

•  In  §2.1,  the  abstract  class  of  constrained  control  problems  that 
we  consider  is  defined. 

•  In  §2.2,  a  list  of  assumptions  about  the  class  of  abstract  prob¬ 
lems  is  given. 

•  In  Theorem  2.1  of  §2.3,  some  of  the  assumptions  listed  in  §2.2 
are  used  to  show  that  optimal  solutions  of  the  abstract  problem 
exist. 

•  In  §2.4,  some  additional  assumptions  of  §2.2  are  used  to  show 
that  Lagrange  multipliers  exist  that  may  be  used  to  enforce  the 
constraint;  also,  first-order  necessary  conditions  are  given. 

•  In  Theorems  2.5  and  2.6  of  §2.4,  the  first-order  necessary  condi¬ 
tions  for  determining  optimal  states  and  controls  are  simplified 
under  additional  assumptions  about  the  control  set. 

•  In  §2.5,  the  optimality  system  from  which  optimal  controls  and 
states  can  be  determined  is  made  more  amenable  to  approxima¬ 
tion  by  simplifying  the  dependence  of  the  objective  functional 
on  the  control. 

2.1.  The  abstract  setting 

We  begin  with  the  definition  of  the  abstract  class  of  nonlinear  control  or  opti¬ 
mization  problems  that  we  study. 

We  introduce  the  spaces  and  control  set  as  follows.  Let  G,  X,  and  Y  be  reflexive 
Banach  spaces  whose  norms  are  denoted  by  ||  •  ||G,  ||  •  ||x  >  and  ||  ■  ||y ,  respectively.  Dual 
spaces  will  be  denoted  by  (•)*.  The  duality  pairing  between  X  and  X*  is  denoted 
by  (■,  -)x  ;  one  similarly  defines  (-,)y  and  {-,)G.  The  subscripts  are  often  omitted 
whenever  there  is  no  chance  for  confusion.  Let  0,  the  control  set,  be  a  closed  convex 
subset  of  G.  Let  Z  be  a  subspace  of  Y  with  a  compact  imbedding.  Note  that  the 
compactness  of  the  imbedding  Z  C  Y  will  play  an  important  role. 

We  assume  that  the  functional  to  be  minimized  takes  the  form 

(2.1)  J(v,z)  =  XT(v)  +  X€{z)  V(i >,z)eXx0, 

where  T  is  a  functional  on  X,  £  a  functional  on  0,  and  A  is  a  given  parameter  which 
is  assumed  to  belong  to  a  compact  interval  A  C  It  + . 

The  constraint  equation  M(v,  z)  =  0  relating  the  state  variable  v  and  the  control 
variable  z  is  defined  as  follow.  Let  AT  be  a  differentiable  mapping  from  X  to  Y,  K  a 
continuous  linear  operator  from  0  to  Y,  and  T  a  continuous  linear  operator  from  Y 
to  X.  For  any  A  6  A,  we  define  the  mapping  M  from  X  x  ©  to  X  by 

(2.2)  M(v,  z)  =  v  +  XTN(v)  +  XTK(z)  V  (v,  z)  €  X  x  0  . 

With  these  definitions  we  now  consider  the  constrained  minimization  problem: 

(2.3)  min  J(v,  z)  subject  to  Af(v,z)  =  0. 

(«,j)exxs 
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In  (2.3),  we  seek  a  global  minimizer  with  respect  to  the  set  {(v,z)  €  X  x  0  : 
M(v,z)  =  0}.  Although,  under  suitable  hypotheses,  we  will  show  that  the  problem 
(2.3)  has  a  solution,  in  practice,  one  can  only  characterize  local  minima,  i.e.,  points 
(ti,  g)  €  X  x  0  such  that  for  some  e  >  0 

J(tt ,  $)  <  J(v,  z)  V  (v,  z)  €  X  x  0  such  that 

M (v,  z)  =  0  and  ||u  —  v||*  <  t . 

Thus,  when  we  consider  algorithms  for  locating  constrained  minima  of  J ,  we  must  be 
content  to  find  local  minima  in  the  sense  of  (2.4). 

After  showing  that  optimal  solutions  exist  and  that  one  is  justified  in  using  the 
Lagrange  multiplier  rule,  we  will  introduce  some  simplifications  in  order  to  render  the 
abstract  problem  (2.3),  or  (2.4),  more  amenable  to  approximation.  The  first  is  to  only 
consider  the  control  set  0  =  G.  The  second  is  to  only  consider  Frechet  differentiable 
functionals  £(•)  such  that  the  Frechet  derivative  £'(g)  =  E~xg,  where  E  is  an  invertible 
linear  operator  from  G*  to  G. 

2.2.  Hypotheses  concerning  the  abstract  problem 

The  first  set  of  hypotheses  will  be  invoked  to  prove  the  existence  of  optimal 
solutions.  It  is  given  by: 

(HI)  infv€x  E(y)  >  —00  ; 

(H2)  there  exist  constants  a,  0  >  0  such  that  £(z)  >  aHzll*5  V  z  E  0 ; 

(H3)  there  exists  a  (v,z)  6  X  x  0  satisfying  M (v,  z)  —  0 ; 

(H4)  if — *  u  in  X  and  g^  g  in  G  where  {(u^"\^n^)}  C  X  x  0,  then 
Af(u("))  _  N(u)  in  Y  and  K(gW)  —  K(g)  tnY  \ 

(H5)  J(-,  -)  is  weakly  lower  semicontinuous  on  X  x  0;  and 

(H6)  if  {(u*”),^"))}  C  X  x  0  is  such  that  {^(u^n^)}  is  a  bounded  set  in  R  and 

M(u(n\g(n))  =  0,  then  {«(")}  is  a  bounded  set  in  X. 

The  second  set  of  assumptions  will  be  used  to  justify  the  use  of  the  Lagrange 
multiplier  rule  and  to  derive  an  optimality  system  from  which  optimal  states  and 
controls  may  be  determined.  The  second  set  is  given  by: 

(H7)  for  each  z  €  0,  v  J(v,z)  and  v  M{y,z)  are  Frechet  differentiable  -, 

(H8)  z  £{z)  is  convex,  i.e., 

£(7*1  +(1  -7)22)  <  7 £(*i)  +  (1  -  7)£(*2)  Vzi,z2  €  0,  V7  €  [0, 1] ; 

and 

(H9)  for  v  €  X,  N'(v)  maps  X  into  Z. 

In  (H9),  N'  denotes  the  Frechet  derivative  of  N. 

A  simplified  optimality  system  may  be  obtained  if  one  invokes  the  additional 
assumption: 

(H10)  0  =  G,  and  the  mapping  z  t— »  £{z)  is  Frechet  differentiable  on  G. 

Hypotheses  (H7)-(H10)  allow  us  to  obtain  a  simplified  optimality  system  for  al¬ 
most  all  values  of  the  parameter  A  €  A.  In  many  cases,  it  is  possible  to  show  that 
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the  same  optimality  system  holds  for  all  values  of  A.  The  following  two  additional 
assumptions  which  will  only  be  invoked  in  case  (1/A)  is  an  eigenvalue  of  -TN'(u) 
each  provides  a  setting  in  which  this  last  result  is  valid: 

(HU)  ifv*  €  X*  satisfies  (/  +  A[W/(u)]*7'*)i;*  =  0  and  K'T'v'  =  0,  then 
v*  =  0 ;  or 

(HI  1 )'  the  mapping  (t>,  z)  >-*  v  +  A TN'(u)v  +XTKz  is  onto  from  X  x  G  to  Y. 

In  order  to  make  the  optimality  system  more  amenable  to  approximation  and 
computation,  we  will  invoke  the  following  additional  assumption: 

(H12)  £'( g )  =  E~lg ,  where-  E  is  an  invertible  linear  operator  from  G *  to  G  and 

g  is  an  optima I  control  for  the  constrained  minimization  problem  (2.4). 


2.3.  Existence  of  an  optimal  solution 


We  first  use  assumptions  (H1)-(H6)  to  establish  that  optimal  solutions  exist. 

THEOREM  2.1.  Assume  that  the  functional  J  and  mapping  M  defined  by  (2.1)  and 
(2.2),  respectively,  satisfy  the  hypotheses  (H1)-(H6).  Then,  there  exists  a  solution  to 
the  minimization  problem  (2.3). 


Proof:  Assumption  (H3)  simply  asserts  that  there  is  at  least  one  element  of  X  x  ©  that 
satisfies  the  constraint.  Thus,  we  may  choose  a  minimizing  sequence  {(«("),  g(n))}  C 
X  x  &  such  that 


lim  J(u(n),0<n>)  = 

n—soo 


inf 

(«i,*)exxe 


J(v,z) 


and 

M(u<n>,s(n>)  =  0. 


By  (HI)  and  (H21,  the  boundedness  of  (,7(t/n*,0(n*)}  implies  the  boundedness  of 
the  sequences  {||</'"^||g}  and  {/■(«("))}.  Then,  by  (H6),  we  deduce  that  is 

bounded.  Thus,  we  may  extract  a  subsequence  {(t/n),  </("))}  such  that  — -  u  in  X 

and  </(")  — -  g  in  G.  Since  0  is  closed  and  convex,  we  have  g  £  0.  Of  course,  u  £  X. 
We  next  show  that  (u,g)  satisfies  the  constraint  equation.  Using  (H4),  we  have  that 

lim  <7W(*<n>),  /)  =  lim  (*(«<">),  T  f)  =  (N(u),  T f)  =  ( TN(u ),  /)  V  /  e  A* 

n— *oq  n— *oo 


and 


lim  (TK(g^),f)  =  lim  (K(gW),T*f)  =  (K{g),T*f)  =  (TK(g)J)  V  /  €  X* 

n  — ♦  oo  n— »oo 

so  that 


0=  lim  (M(u<">, *<">),  /}  =  («  +  A  7W(u)  +  A TK (<?),/)  V  /  G  X*  , 

n— *oo 

i.e.,  M(u,g)  =  0.  Finally,  we  use  (H5),  the  weak  lower  semicontinuity  of  J(-,  •),  to 
conclude  that  ( u,g )  is  indeed  a  minimizer  in  X  x©  satisfying  the  constraint  M(u,g)  = 
0.  □ 


2.4.  Existence  of  Lagrange  multipliers 

We  now  wish  to  use  the  additional  assumptions  (H7)-(H9)  to  show  that  the  La¬ 
grange  multiplier  rule  may  be  used  to  turn  the  constrained  minimization  problem 
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(2.3)  into  an  unconstrained  one.  Actually,  the  Lagrange  multiplier  rule  will  only  en¬ 
able  us  to  find  local  minima  in  the  sense  of  (2.4).  We  first  quote  the  following  abstract 
Lagrange  multiplier  rule  whose  proof  can  be  found  in  [21]. 

THEOREM  2.2.  Let  X\  and  X 2  be  two  Banach  spaces  and  0  an  arbitrary  set.  Suppose 
J  is  a  functional  on  X\  x  0  and  M  a  mapping  from  X\  x  0  to  X2.  Assume  that 
(u,  9)  €  X\  x  0  is  a  solution  to  the  following  constrained  minimization  problem: 

M(u,g)  =  0  and  there  exists  an  e  >  0  such  that  J(u,g )  <  J(v,z) 

*  for  all  ( v ,  z)  €  Xi  x  0  such  that  ||u  —  v||x,  <  f  and  M(v,  z)  =  0 . 

Let  U  be  an  open  neighborhood  of  u  in  X\ .  Assume  further  that  the  following  conditions 
are  satisfied: 

(2.6)  for  each  z  €  0,  v  *-*  J{y ,  z)  and  v  >-*  M(y,  z)  are  Frichet- differentiable  at 
v  =  u ; 

(2.7)  for  any  v  €  U,  z\,z2  €  0,  and  y  6  [0, 1],  there  exists  a  z7  =  z7(t>, zi,zz) 
such  that 

M(v,  z7)  =  yM(v,  zi )  +  ( 1  -  y)M (t>,  z2) 

and 

<  yJ(v,  z,)  +  (1  -  y)J(v,  z2) ; 

and 

(2.8)  Range  (Mu(u, g))  is  closed  with  a  finite  codimension, 

where  Mu(u,g)  denotes  the  Frechet  derivative  of  M  with  respect  to  u.  Then,  there 
exists  a  k  £  R  and  a  /i  €  X2  that  are  not  both  equal  to  zero  such  that 

*  (Ju(u,  g),  v)  -  (/i,  M„(«,  g)v)  =  0  V»€^i 

and 

min  C(u,  z,  n,  *)  =  C(u,  g,  n,  *) , 

where  C(u,g,n,k)  =  kj(u,g )  —  ( n,M(u,g ))  is  the  Lagrangian  for  the  constrained 
minimization  problem  (2.5)  and  where  Ju{n,g)  denotes  the  Frechet  derivative  of  J 
with  respect  to  u.  Moreover,  if 

(2.9)  the  algebraic  sum  Mu(u,g)X\  +  Af(u,0)  contains  0  €  X2  us  an  interior 
point, 

then  we  may  choose  k  =  1,  i.e.,  there  exists  a  n  €  X2  such  that 
{Ju{u,g),v)  -  (p,  Mu(u,g)v)  =  0  Vt)6Xi 

and 

min£(u,  z,ft,l)=  C(u,g,n,  1) . 

Proof:  See  [21].  0 

Next,  we  apply  Theorem  2.2  to  the  optimization  problem  (2.4).  In  doing  so,  we 
will  need  the  following  result. 

LEMMA  2.3.  Let  the  spaces  X,  Y ,  and  Z  and  operators  T  and  N  be  defined  as  in  §2.1. 
For  v  €  X,  assume  that  N'(v)  maps  X  into  Z.  Then,  TN'{y)  is  a  compact  operator 
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from  X  to  X,  and  therefore  <t(—TN'(v)),  the  spectrum  of  the  operator  (—TN'(v)),  is 
at  most  countable  with  zero  being  the  only  possible  limit  point. 

Proof:  Since  Z  Y,  we  see  that  N'(v)  is  a  compact  linear  operator  from  X  to 
Y.  Also,  T  is  a  bounded  linear  operator  from  V  to  A-  so  that  TN'(v)  is  a  compact 
operator  from  X  to  X.  Hence,  ff(—TN'(v))  is  at  most  countable  and  consists  only  of 
0  and  the  eigenvalues  of  (- TN'(v )) .  □ 

Note  that  in  the  following  result,  the  existence  of  at  least  one  pair  (u,  g)  satisfying 
(2.4)  is  guaranteed  by  Theorem  2.1. 

THEOREM  2.4.  Let  A  €  A  be  given.  Assume  that  assumptions  (H1)-(H9)  hold.  Let 
(u,g)  €  X  x  Q  be  an  optimal  solution  satisfying  (2.4).  Then,  there  exists  a  k  G  R 
and  a  ft  €  X*  that  are  not  both  equal  to  zero  such  that 

(2.10)  k(Ju(u,g),w)  -  (ft,Mu(u,g)w)  =  0  VtogA 

and 

(2.11)  nun  C(u,  z,  ft,  k)  =  £(ti,  g,  ft,k). 

*6“ 

Furthermore,  i/(l/A)  ^  <r(-TN'(u)),  we  may  choose  1:1,  i.e.,  for  almost  all  X, 
there  exists  a  ft  G  X*  such  that 

(2.12)  g),  w)  -  {ft,  M„( A,  u,  g)  ■  w)  =  0  V  w  6  X 

and 

(2.13)  min£(ti,  z,ft,l)  =  £(u,  g,  ft,  1). 


Proof:  Let  A  €  A  be  given.  To  show  the  existence  of  k  and  ft  such  that  (2.10)  and 
(2.11)  are  valid,  we  only  need  to  verify  that  the  hypotheses  (2.6)-(2.8)  of  Theorem  2.2 
hold  with  X\  =  Xi  —  X,  since  in  this  case  (2.5)  reduces  to  (2.4).  Obviously,  (2.6)  is 
merely  a  restatement  of  (H7).  Since  8  is  convex  and  since  the  mappings  T  and  K  are 
linear,  we  have  that  if  z-,  =  yz\  +  (1  —  7)22.  then 

Af(v,Zy)  =  v  +  XTN(v)  +  A  TKz^ 

=  y(v  +  XTN{v))  +  (1  -  7)(v  +  A TN{ v))  +  yX {TKz,  +  (1  -  y)TKz2)) 

=  yM(v,zi)  +  (l-y)M(v,z2). 

Moreover,  (H8)  implies  that 

J(v,Zy)  =  A^(v)  +  XS(z-f)  =  XT{v)  +  A  5(721  +  (1-7)22) 

<  A/-(u)  + A  (75(21) +  (1  -7)5(22))  =7  J(v,2i)  +  (1  -y)J{v,z2). 

Thus,  (2.7)  holds.  The  operator  Mu(u,g)  from  X  to  X  is  defined  by 

Mu(v,  g)  ■  w  =  w  +  A7W'(u)  •  u;  VttgX 


or  simply, 


Mu(u,g)  =  I  +  XTN'{u)  . 


From  (H9)  and  Lemma  2.3,  we  have  that  TN'( u)  is  a  compact  operator  from  X  to  X. 
As  a  result,  Mu(u,g)  =  /  +  A TN'(u)  is  a  Fredholm  operator  so  that  it  has  a  closed 
range  with  a  finite  codimension,  i.e.,  (2.8)  holds.  Thus,  by  Theorem  2.2,  there  exists 
a  k  €  R  and  a  n  €  X*  which  are  not  both  equal  to  zero  such  that  (2.10)  and  (2.11) 
hold. 

To  show  the  existence  of  a  n  such  that  (2.12)  and  (2.13)  are  valid,  we  only  need 
to  verify  that  the  additional  hypothesis  (2.9)  of  Theorem  2.2  holds.  In  fact,  if,  in 
addition  (1/A)  g  <t(—TN'(u)),  then  it  follows  that  X  =  Range(7  +  A TN'(u))  = 
Range(A/„(u,  (?))  so  that  Range(Mu(u,y))  contains  0  €  A  as  an  interior  point,  i.e., 
(2.9)  holds.  Hence,  by  Theorem  2.2  and  Lemma  2.3,  we  conclude  that  for  almost  all 
A,  there  exists  a  /i  €  X *  such  that  (2.12)  and  (2.13)  hold.  □ 

So  far  0  has  only  been  assumed  to  be  a  closed  and  convex  subset  of  G.  No 
smoothness  condition  on  the  control  variable  g  has  been  assumed  in  the  functional 
or  in  the  constraint.  Thus,  the  necessary  condition  of  optimality  with  respect  to 
variations  in  the  control  variable  is  expressed  in  the  cumbersome  relation  (2.11).  We 
now  turn  to  the  case  where  0  contains  a  neighborhood  of  g,  where  (u,g)  is  an  optimal 
solution.  In  particular,  we  assume  that  0  =  G.  In  this  case,  (2.11)  can  be  given  a 
more  concrete  structure. 

THEOREM  2.5.  Let  A  e  A  be  given.  Assume  that  assumptions  (H1)-(H10)  hold.  Let 
(ti,  g)  €  X  x  G  be  a  solution  of  the  problem  (2.4).  Then,  there  exists  a  k  G  R  and  a 
H  €  X*  that  are  not  both  equal  to  zero  such  that 

(2.14)  •  k  {Ju(u,  g),w)-  (p,  (I  +  XT  N'(u))w)  =  0  VuigA 

and 

(2.15)  k(€'(g),  z)  -  (/i,  TKz)  =0  V  x  €  G . 

Furthermore,  i/(l/A)  g  <r(— TN'(u)),  we  may  choose  k  =  1,  i.e.,  there  exists  a  n  6  X" 
such  that 

(2.16)  (Ju(u,g),w)  —  (fz,  (/  +  \TN'(u))w)  =  0  V»eA 
and 

(2.17)  (6'{g),z)  —  (fx,TKz)  =  0  Vz<=G. 
hold. 

Proof:  Since  the  hypotheses  imply  that  J(y,z)  is  Frechet  differentiable  with  respect 
to  z,  (2.14)-(2.17)  follow  easily  from  Theorem  2.4.  □ 

Remark.  If  Jb  =  0,  then  there  exists  a  n  ^  0  such  that 

-(/i,  Mu(u,  g)w)  =  0  Vine  A 

so  that  the  optimality  system  necessarily  has  infinitely  many  solutions.  In  fact,  for  any 
C  €  R ,  {Cp)  is  a  solution  whenever  /i  is  a  solution.  This  creates  both  theoretical  and 
numerical  difficulties.  Thus,  it  is  of  great  interest  to  try  to  eliminate  this  situation. 
Fortunately,  Lemma  2.3  and  Theorem  2.4  tell  us  that  we  may  set  k  =  l  ^  0  for  almost 
all  values  of  (1/A),  i.e.,  except  for  the  at  most  countable  set  of  values  in  <r(— TN’(u)). 
□ 
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If  the  control  g  enters  the  constraint  in  a  favorable  manner,  then  we  may  take 
Jb  =  1  even  when  (1/A)  €  <r(— TJV'(u)).  Specifically,  we  invoke  one  of  the  assumptions 
(HU)  and  (Hll)'.  We  then  have  the  following  result. 

THEOREM  2.6.  Assume  that  the  hypotheses  of  Theorem  2.5  hold.  Assume  that  if 
(1/A)  €  <r(—TN'(u)),  then  either  (Hll)  or  (Hll)'  holds.  Then,  for  all  A  €  A,  there 
exists  a  p  €  X*  such  that  (2.16)  and  (2.17)  hold. 

Proof:  Because  of  Theorem  2.5,  we  only  need  to  examine  the  case  (1/A)  €  <r(— T  N'(u)) 
and  show  that  the  algebraic  sum  Mu(u,g)X-rM(u,  G)  =  X.  If  (Hll)'  holds,  the  result 
is  a  direct  application  of  Theorem  2.2. 

If  (Hll)  holds,  let  (1/A)  be  a  nonzero  eigenvalue  of  (-TN'(u)).  Then,  A  is 
also  an  eigenvalue  of  (—JV'(u)0T*)  with  a  finite  dimensional  eigenspace  having  the 
corresponding  eigenfunctions  (ti*  }"Lj  C  X*  as  a  basis.  We  claim  that  {K’T* v* ,  C 
G*  is  a  linearly  independent  set.  To  see  this,  we  assume  j  a,  K’T*v *  =  0  with 
a,  e  R;  this  expression  can  be  rewritten  as  K’T*  (E,™  i  a,  t>* )  =  0.  Because  each  t>* 
is  an  eigenvector,  we  have  (/  +  A7V'(u)*7”)  EcLi  a«  vi  —  0-  Thus,  the  assumption 
(Hll)  implies  that  Ei*LiQ(«t’»*  ==  0-  Since  is  an  eigenbasis,  and  is  therefore 

a  linearly  independent  set,  we  have  each  ati  =  0.  This  shows  that  {K'T’ u*  JJ’lj  is 
linearly  independent  set  in  G*.  Hence,  we  may  choose  an  orthonormal  dual  basis 
{*.}£i  C  G  such  that  (zi,K*T*vj)  =  Sij. 

Now,  let  w  €  X  be  given.  We  choose  z  =  x YlT=i(w’vi)z>-  Then  (u>,t>?)  — 
A (TKz,  t,;)  =  (w,  V »J)  -  A(z,  K*T-V))  =  (w,  t/J)  -  E^iK «?)«y  =  0  for  j  =  1, . . . ,  m. 
Thus,  by  Fredholm  alternatives,  there  exists  a  unique  v  £  X  that  satisfies  (/  + 
XTN'(u))v  =  w  -  A  TKz,  or,  (I  +  XTN'(u))v  +  \TKz  =  w ;  thus,  we  have  shown 
that  Mu(u,g)X  +  M(u,  G)  =  X.  Hence,  by  Theorem  2.2,  there  exists  a/i£A*  such 
that  (2.16)  and  (2.17)  hold.  □ 

2.5.  The  optimality  system 

Under  the  assumptions  of  Theorem  2.6,  an  optimal  state  u  G  X,  an  optimal  con¬ 
trol  g  €  G,  and  the  corresponding  Lagrange  multiplier  /j  €  X*  satisfy  the  optimality 
system  of  equations  formed  by  (2.2),  (2.16),  and  (2.17).  From  (2.1)  we  have  that 
Ju  =  XT'  and  J}  —  XS',  where  T'  denotes  the  obvious  Frechet  derivative.  Then, 
(2.16)-(2.17)  may  be  rewritten  in  the  form 

(2.18)  n  +  X  [N'(u)YT'n  -  X  T\ u)  =  0  in  X* 
and 

(2.19)  Z\g)  -  K*T* fi  =  0  in  G* . 

For  purposes  of  numerical  approximations,  it  turns  out  to  be  convenient  to  make  the 
change  of  variable  £  =  T*p.  Then,  the  optimality  system  (2.2),  (2.18),  and  (2.19)  for 
u  €  X ,  g  €  G,  and  £  €  Y*  takes  the  form 

(2.20)  ti  +  A7W(ti)  +  XT  Kg  =0  in  X  , 

£  +  A T*[N'(u)]*£  —  XT*T'(u)  =  0  in  Y*  , 


€\g)  -  K*Z  —  0  in  G* . 


(2.21) 

and 

(2.22) 
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It  will  also  be  convenient  to  invoke  an  additional  simplifying  assumption  concern¬ 
ing  the  dependence  of  the  objective  functional  on  the  control.  Specifically,  we  assume 
that  (H12)  holds.  Then,  (2.20)-(2.22)  can  be  rewritten  as 


(2.23) 

u+\TN(u)  +  ATKg  =  0  in  X  , 

(2.24) 

(  +  AT*[N'(u)]'(-  ir^(u)  =  0  in7* 

and 

(2.25) 

g  -  EK’(  =  0  inG. 

Remark.  Note  that  the  optimality  systems,  e.g.,  (2.23)-(2.25),  are  linear  in  the  adjoint 
variable  £.  Also,  note  that  the  control  g  may  be  eliminated  from  the  optimality  system 
(2.23)-(2.25).  Indeed,  the  substitution  of  (2.25)  into  (2.23)  yields 

(2.26)  u  +  A  TN(u)  +  A  TKEK't  =  0  in  X  . 

Thus,  (2.24)  and  (2.26)  determine  the  optimal  state  u  and  adjoint  state  £;  subse¬ 
quently,  (2.25)  may  be  used  to  determine  the  optimal  control  g  from  £.  This  observa¬ 
tion  serves  to  emphasize  the  important,  direct  role  that  the  adjoint  state  plays  in  the 
determination  of  the  optimal  control.  □ 

Remark.  Given  a  £  G  Y*,  it  is  not  always  possible  to  evaluate  g  exactly  from  (2.25). 
For  example,  the  application  of  the  operator  E  may  involve  the  solution  of  a  partial 
differential  equation.  Thus,  although  it  is  often  convenient  to  devise  algorithms  for 
the  approximation  of  optimal  control  and  states  based  on  the  simplified  optimality 
system  (2.24)  and  (2.26),  in  some  other  cases  it  is  best  to  deal  with  the  full  form 
(2.23)-(2.25).  Thus,  when  we  consider  approximations  of  optimal  controls  and  states, 
we  will  deal  with  the  latter.  □ 

Remark.  In  many  applications  we  have  that  X*  =  Y.  Since  these  spaces  are  assumed 
to  be  reflexive,  we  also  have  that  Y *  =  X.  In  this  case,  we  have  that  both  it  and  £ 
belong  to  X.  □ 


3.  Finite  dimensional  approximations  of  the  abstract  problem 

In  this  section  we  define  and  analyze  algorithms  for  the  finite  dimensional  approx¬ 
imation  of  solutions  of  the  optimality  system  (2.23)-(2.25);  an  outline  of  the  definitions 
and  results  of  this  section  is  as  follows. 

•  In  §3.1,  we  define  the  finite  dimensional  approximate  problems 
that  we  consider. 

•  In  §3.2,  a  list  of  assumptions  about  the  approximate  problems 
is  given. 

•  In  §3.3,  we  quote  a  result  of  [6]  that  we  will  use  to  analyze 
approximations  obtained  as  solutions  of  the  approximate  prob¬ 
lems  defined  in  §3. 1-3.2. 

•  In  §3.4,  we  provide  error  estimates  for  the  approximation  of 
solutions  of  the  optimality  system  (2. 23)- (2. 25) 
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3.1.  Formulation  of  finite  dimensional  approximate  problems 

A  finite  dimensional  discretization  of  the  optimality  system  (2.23)-(2.25)  is  de¬ 
fined  as  follows.  First,  one  chooses  families  of  finite  dimensional  subspaces  Xh  C  X , 
(y»)h  £  y^  an(j  Qh  q  These  families  are  parameterized  by  a  parameter  h  that 
tends  to  zero.  (For  example,  this  parameter  can  be  chosen  to  be  some  measure  of  the 
grid  size  in  a  subdivision  of  Cl  into  finite  elements.)  Next,  we  define  approximate  oper¬ 
ators  T*  :  Y  -*  Xh,  Eh  :  G*  —  GA,  and  (T‘)h  :  X *  —  (y*)A.  Of  course,  one  views 
T\  EA,  and  (T*)A  as  approximations  to  the  operators  T,  E,  a.  J  7” ,  respectively. 
Note  that  (T*)A  is  not  necessarily  the  same  as  (TA)*.  The  former  is  a  discretization 
of  an  adjoint  operator  while  the  later  is  the  adjoint  of  a  discrete  operator. 

Once  the  approximating  subspaces  and  operators  have  been  chosen,  an  approxi¬ 
mate  problem  is  defined  as  follows.  We  seek  uh  £  Xh,  gh  £  GA,  and  £  (Y’)h  such 
that 

(3.1)  uh  +  XThN{uh)  +  XThKgh  =0  in  Xh  , 

(3.2)  £h  +  X(T*)li[N'(ull)]',th  —  A(T*)h^’'(ufc)  =  0  in  (y*)A  , 
and 

(3.3)  gh  -  EhKmth  =  0  in  Gh  . 


3.2.  Hypotheses  concerning  the  abstract  problem  and  the  approximate 
problem 

We  make  the  following  hypotheses  concerning  the  approximate  operators  Th , 
(T*)A,  and  Eh  \ 

(H13)  lim||(T-Th)y|U  =0  Vygy, 

rl-*0 

(H14)  lim  \\(T*  —  (T*)A)i;||y  =  0  V  v  £  X*  , 

h— *0 

and 

(H 15)  lim  ||(£  -  £A)s||G  =0  V  s  €  G* . 

/i— *0 

We  also  need  the  following  additional  hypotheses  on  the  operators  appearing  in 
the  definition  of  the  abstract  problem  (2.4): 

(H16)  N  £C3(X,Y)  and  F  £C3(X-,R)  ] 

(H17)  N" ,  N'",  T" ,  and  T'"  art  locally  bounded,  i.e.,  they  map  bounded  sets  to 
bounded  sets; 

(H18)  for  v  £  X,  in  addition  to  (H9),  i.e.,  N'(v)  £  C(X\Z)  where  Z  Y , 
we  have  that  [N'(«)]*  €  C{Ym\Z)  where  Z  <-*<—>  X* ,  that  for  t]  £  Y* , 
[AT(v)]‘  •  r}  £  C(Y*;Z),  and  that  for  w£X,  T"{v)  •  w  £  C{X ;  Z);  and 
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(H19)  K  maps  G  into  Z. 

Here,  (■)"  and  (■)"'  denote  second  and  third  FYechet  derivatives,  respectively. 

3.3.  Quotation  of  results  concerning  the  approximation  of  a  class  of  non¬ 
linear  problems 

The  error  estimate  to  be  derived  in  Section  3.4  makes  use  of  results  of  [6]  and  [10] 
(see  also  [13])  concerning  the  approximation  of  a  class  of  nonlinear  problems.  These 
results  imply  that,  under  certain  hypotheses,  the  error  of  approximation  of  solutions 
of  certain  nonlinear  problems  is  basically  the  same  as  the  error  of  approximation  of 
solutions  of  related  linear  problems.  Here,  for  the  sake  of  completeness,  we  will  state 
the  relevant  results,  specialized  to  our  needs. 

The  nonlinear  problems  considered  in  [6],  [10],  and  [13]  are  of  the  following  type. 
For  given  A  €  A,  we  seek  tfi  G  X  such  that 

(3.4)  H(\,rP)  =  tP  +  TG(\,v)  =  0, 

where  T  G  C(y-,X),  G  is  a  C2  mapping  from  A  x  X  into  X  and  y  are  Banach 
spaces,  and  A  is  a  compact  interval  of  R .  We  say  that  {(A,  t/j(A))  :  A  G  A}  is  a  branch 
of  solutions  of  (3.4)  if  A  — ►  ^>(A)  is  a  continuous  function  from  A  into  X  such  that 
7f(A,  V’(A))  =  0.  The  branch  is  called  a  regular  branch  if  we  also  have  that  W^(A,  t/»(A)) 
is  au  isomorphism  from  X  into  X  for  all  A  G  A.  Here,  •)  denotes  the  Frechet 

derivative  of  W(-,  -)  with  respect  to  the  second  argument.  We  assume  that  there  exists 
another  Banach  space  Z,  contained  in  y,  with  continuous  imbedding,  such  that 

(3.5)  g^X  ,J>)eC(X-,Z)  VAgA  and  ip  €  X , 

where  Gy(-,  •)  denotes  the  Frechet  derivative  of  {/(•,  )  with  respect  to  the  second  ar¬ 
gument. 

Approximations  are  defined  by  introducing  a  subspace  Xh  C  X  and  an  approxi¬ 
mating  operator  Th  G  C(y,Xh).  Then,  given  A  G  A,  we  seek  iph  G  Xh  such  that 

(3.6)  Hh(\,1>h)  =  il>h+ThG(X,tl>h)  =  0. 

Concerning  the  operator  Th ,  we  assume  the  approximation  properties 

(3.7)  lim  ||(Ta  -  7>||*  =0  V  w  G  y 

h—*0 

and 

(3.8)  lim||(Tfc-T)||£(^)  =  0. 

Note  that  whenever  the  imbedding  Z  C  ^  is  compact,  (3.8)  follows  from  (3.7)  and, 
moreover,  (3.5)  implies  that  the  operator  TG${ X,rp)  G  C(X,X)  is  compact. 

We  can  now  state  the  result  of  [6]  or  [10]  that  will  be  used  in  the  sequel.  In  the 
statement  of  the  theorem,  D2G  represents  any  and  all  second  Frechet  derivatives  of  Q. 

THEOREM  3.1.  Let  X  and  y  be  Banach  spaces  and  A  a  compact  subset  of  R.  Assume 
that  G  is  a  C2  mapping  from  Ax  X  into  y  and  that  D2G  is  bounded  on  all  bounded 
sets  of  A  x  X.  Assume  that  (3.5),  (3.7),  and  (3  8)  hold  and  that  {(A,  V’M);  A  G  A)  is  a 
branch  of  regular  solutions  of  (3.4).  Then,  there  exists  a  neighborhood  O  of  the  origin 
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in  X  and,  for  k  <  ho  small  enough,  a  unique  C2  function  A  — *  tph(X)  G  Xh  such  that 
{(A,  0fc(A));  A  €  A}  is  a  branch  of  regular  solutions  of  (3.6)  and  -  rp( A)  G  O  for 
all  A  €  A.  Moreover,  there  exists  a  constant  C  >  0,  independent  of  h  and  A,  such  that 

(3.9)  ||V.k(A)-V'(A)|U<C||(Tfc-T)C(A,V-(A))|U  VAgA.  □ 

3.4.  Error  estimates  for  the  approximation  of  solutions  of  the  optimality 
system 

We  now  apply  the  result  of  Theorem  3.1  to  study  the  approximation  of  solutions 
of  the  optimality  system.  Set  X  =  X  x  G  x  Y *,  y  =  Y  x  X *,  Z  =  Z  x  Z,  and 
Xh  =  Xh  xCkx  ( Y*)h .  (Recall  that  Z  was  introduced  in  (H18).)  By  the  hypotheses 
on  Z  and  Z,  we  have  that  Z  is  compactly  imbedded  into  y.  Let  T  G  £(>’;A’)  be 
defined  in  the  following  manner:  T(r,  f)  =  (u,jj.O  for  (r,f)  G  y  and  {u,g,Q  G  X  if 
and  only  if 

(3.10)  u  +  Tf  =  0, 

(3.11)  £  +  r*f  =  o, 

and 

(3.12)  g-EK’i=  0. 

Similarly,  the  operator  Tfc  G  C(y;Xh)  is  defined  as  follows:  Th(f,f)  =  (uh,ph,|h) 
for  (r,  f )  G  >  and  (uh,  gh,th)  G  if  and  only  if 


(3.13) 

uh  +  Thf  =  0 , 

(3.14) 

ih  -f  (T*)hf  =  0, 

and 

(3.15) 

gh-EhK*ih  =  0. 

The  nonlinear  mapping  Q  :  A  x  X  — *■  y  is  defined  as  follows:  f?(A,(u ,g,£)) 
for  A  G  A,  (ti,  g,£)  €  X,  and  (f,  f)  G  y  if  and  only  if 

(3.16) 

r  =  \  N{v)  +  XKg 

and 

(3.17)  f 

=  A  [jV'(u)]*£  -  A  T'{u) . 

It  is  evident  that  the  optimality  system  (2.23)-(2.25)  and  its  finite  dimensional  coun¬ 
terpart  (3-l)-(3.3)  can  be  written  as 

( u,g,O  +  Tg(X,(u,g,O)  =  0 
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and 


(n\gh,e)  +  Thg{X,(n\gh,ih))  =  0, 

respectively,  i.e.,  with  rp  =  (t i,g,£)  and  iph  =  (t tk,gk,il>k),  in  the  form  of  (3.4)  and 
(3.6),  respectively. 

Now  we  examine  the  approximation  properties  of  Tk. 

LEMMA  3.2.  Let  the  operatorsT  and  Tk  be  defined  by  (3.10)-(3.12)  and  (3. 13)-(3. 15), 
respectively.  Assume  that  the  hypotheses  (H13)-(H15)  hold.  Then, 

(3.18)  lim||(T-TA)(r,r)(U  =  0  V(r  ,r)G^. 

A— *0 


Proof:  Let  (u,ij,0  =  T(r,r),  i.e.,  («,£,£)  satisfies  (3.10)-(312).  Let  (uh,gh,(-h)  = 
TA(r,r)  ,  i.e.,  (uk,gk,^k)  satisfies  (3.13)-(3.15).  Subtracting  the  corresponding  equa¬ 
tions  yields  that 

||u-uA|U=||(T-rA)r||x, 

IK-«kl|y  =  ll(T*-(nk)r||v. 

and 

119  ~  9h\\a  =  ||(fi  -  Eh)K'th  +  EK*(t  -  f>)||G 

<  || (E  -  Ek)K^k\\G  +  \\EK'\\c(y.;g)  ||«  -  ik)\ [a  ■ 

Thus,  for  some  constant  C  >  0, 

||(T-Th)(r,r)|U 

<  C  {||(T  -  TA)r||x  +  ||(T*  -  (T*)A)r||y.  +  ||(£  -  £7fc) ||g }  • 


Then,  the  result  of  the  proposition  follows  from  (H13)-(H15).  □ 

Next,  we  examine  the  derivative  of  the  mapping  Q. 

LEMMA  3.3.  Let  the  mapping  Q:\xX  —*  y  be  defined  by  (3.16)-(3.17).  Assume  that 
the  hypotheses  (H9),  (H16),  and  (H18)-(H19)  hold.  Then,  for  every  A  g  A  and  every 
(«.Sf,OeT,  the  operator  £(«,,,<)  (M«>  9.0)  €  C(X;Z). 

Proof:  A  simple  calculation  shows  that  ^(«,j,e)(A ,(«,?,£))  £  £(T;^0  is  given  by 


£(«,#.<)  (A>(«>  9.0)  («,9,0 


N'(u)u+Kg 

[N"(u)  u]*  •£  +  [JV'(u)J*  •  i -  /•"(«) 


-.)■ 


Then,  the  result  follows  from  (H9)  and  (H18)-(H19).  □ 

A  solution  g(\),£(X))  of  the  optimality  system  (2.23)-(2.25)  is  called  regular 

if  the  system  (for  the  unknowns  («,</,£)) 

(3.19)  u  +  A  TN'(u)u  +  A  TK g  =  x , 


(3.20)  £  +  AT*[N"(u)]*u  ^  +  Ar[lV'(u)]*^-  \TT"(u)u  =  y 
and 

(3.21)  g  —  EK’i  =  z  . 
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is  uniquely  solvable  for  any  ( x,z,y )  6  X  =  X  x  G  x  Y* .  (Note  that  the  linear 
operator  appearing  on  the  left  hand  side  of  (3.19)-(3.21)  is  obtained  by  linearizing  the 
optimality  system  (2.23)-(2.25)  about  (u,j,£).) 

In  the  following  theorem,  we  will  assume  that  the  solution  (u(A),  g( A),£(A))  of  the 
optimality  system  (2.23)-(2.25)  that  we  are  trying  to  approximate  is  a  regular  solution. 
The  assumptions  we  have  made,  in  particular  (H9),  (H18)-(H19),  are  sufficient  to 
guarantee  that  for  almost  all  values  of  A,  this  is  indeed  the  case. 

LEMMA  3.4.  Assume  the  hypotheses  of  Lemma  3.3.  Then,  for  almost  all  A,  solutions 
(ti(A),<;(A),£(A))  of  the  optimality  system  (2.23)-(2.25)  ore  regular. 

Proof:  The  system  (3.19)-(3.21)  is  equivalent  to 

(3.22)  (l  +  \TS(u,g,t))(u,g,i)  =  (x,z,y), 
where  the  linear  operator  S(u,g,£) :  X  — *  y  is  defined  by 

S(u,0,£)  •(«,$,£)  =  |5(u,,,e)(A,(u,j,0)  (u,g,£) 

_  (  N'(u)  ■  u  +  Kg  \ 

-  V  [W»  uj*  •{  +  [N'{u)Y  ■  £  -  ?"{u)  u)  ■ 

Now,  T  €  C(y,X),  so  that,  by  Lemma 3.3,  (/  +  \TS(u,g,£))  is  a  compact  perturba¬ 
tion  of  the  identity  operator  from  X  to  X.  Thus,  for  almost  all  A,  (3.22),  or  equivalently 
(3.19)-(3.21),  is  uniquely  solvable,  i.e.,  for  almost  all  A,  the  solution  (u(A),  ff(A),£(A)) 
of  the  optimality  system  (2.23)-(2.25)  is  regular.  □ 

Using  Theorem  3.1,  we  can  now  provide  an  error  estimate  for  appr  .mations  of 
solutions  of  the  abstract  problem. 

THEOREM  3.5.  Let  (u(A),p(A),£(A))  e  X,  for  A  e  A ,  be  a  branch  of  regular  solutions 
of  the  optimality  system  (2.23)-(2.25).  Assume  that  the  hypotheses  (H13)-(H19)  hold. 
Then,  there  exists  a  6  >  0  and  an  ho  >  0  such  that  for  h  <  h0,  the  discrete  optimality 
system  (3.1)-(3.3)  has  a  unique  solution  (uh(A), jrh(A),£h(A))  satisfying 

ll(«(A),*(A),*(A))  -  (u*(A),/(AUk(A))||*  <6. 

Moreover, 

(3.23)  Um||(t.(A),?(A),((A))  -  (u*(A),/(A),(*(A))||*  =  0 

uniformly  in  A  €  A  and  there  exists  a  constant  C,  independent  of  h  and  A,  such  that 

lim  ||(«(A),0(A),C(A))  -  (uft(A),/(A),^(A))|U 

(3.24)  <  CA{|J(Th  -  T){N(u( A))  +  K<,(A))|)X  + 1 \(Eh  -  E)K'Z( A)||G 

+ 1|  ((T*)fc  -  r )  ([^'(u(A))]*e  -  ^'(u( A)))  ||y. }  . 


Proof:  Assumptions  (H16)  and  (H17)  ensure  that  Q  €  C2(X,y)  and  D7Q  maps 
bounded  sets  of  A  x  X  into  bounded  sets  of  y.  By  Lemma  3.3,  assumptions  (H 18) 
and  (H19)  imply  that  (3.5)  holds.  By  Lemma  3.2,  assumptions  (H13)-(H15)  imply 


14 


I 


that  (3.7)  holds.  Then,  since  Z  is  compactly  imbedded  into  y,  (3.7)  implies  that  (3.8) 
holds.  Thus,  all  the  hypotheses  of  Theorem  3.1  are  verified.  Then,  a  direct  application 
of  Theorem  3.1  yields  (3.23)  and  (3.24)  follows  from  (3.9).  □ 

It  is  easily  seen  that  (3.23)  and  (3.24)  are  equivalent  to: 

lim{||u(A)  -  uk(X)\\x  +  ||«(A)  -  9h(X)\\a  +  IMA)  -  *fc(A)||y.}  =  0 

uniformly  in  A  €  A  and  that  there  exists  a  constant  C,  independent  of  h  and  A,  such 
that 

HA)  -  «k(A)||*  +  IMA)  ~ /(A)||c  +  IMA)  -  *k(A)||y. 

<  CA{||(TA  -  T)(N{n( A))  +  tf*(A))||x  +  ||(£k  -  £)***(A)||g 
+  ||((r  )k  -  r)([tf'(«(A))N(A)  -  r(u(A)))||y.} . 

If,  in  (3.9),  the  operator  T  is  invertible,  we  have,  using  (3.4),  that 
IIV'fc(A)  -  tHA)||*  <  C\\(7hT-1  -  7MA)|U  V  A  €  A  . 

Thus,  if  the  operator  T  from  Y  to  X  is  invertible,  we  have  that  (3.24)  is  equivalent  to 
||«(A)  -  uk(A)||*  +  IMA)  -  ffh(A)||c  +  IMA)  -  £k(A)||y* 

(3.25)  <  cjlKT^r-1  -  /)«(A)||x  +  ||(^£r>  -  /)s(A)||g 

+„((n*(n_,-/)«A)iiy}. 


4.  Applications 

We  now  apply  the  framework  and  analyses  developed  in  §2  and  §3  to  some  concrete 
problems,  all  of  which  feature  constraints  on  admissible  states  and  controls  that  take 
the  form  of  a  system  of  nonlinear  partial  differential  equations.  In  each  application, 
we  use  a  different  control  mechanism  so  that  the  discussion  provided  in  this  section 
illustrates  the  treatment  of  a  variety  of  such  mechanisms.  However,  one  could  use  any 
of  the  control  mechanisms  discussed  in  any  of  the  applications  in  any  other  application, 
or  in  fact,  use  any  combination  of  such  mechanisms. 

Before  examining  any  specific  application,  we  establish  some  notation.  Further 
notation  will  be  established  as  needed  when  the  individual  applications  are  considered. 

Throughout,  C  will  denote  a  positive  constant  whose  meaning  and  value  changes 
with  context.  Also,  H'( V)  for  s  6  R  denotes  the  standard  real  Sobolev  space  of  order 
s  with  respect  to  the  set  V,  where  V  could  either  be  a  bounded  domain  Q  G  Rd, 
d  =  2,3,  or  part  of  the  boundary  T  of  such  a  domain.  Of  particular  interest  are  the 
spaces  H°(V)  =  L2( V), 

H\V)={<t>€V\V)  |  ^ei!( V)  for  j  =  1, . .  .,d| 

and 

H’(D)  ={*£*’(!>)  | 
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Abo  of  interest  is  the  subspace 

//,}(2>)  =  {*€//*(*>)  |  ^  =  0  on£2>}, 
where  dV  denotes  the  boundary  of  V. 

Dual  spaces  will  be  denoted  by  (-)*.  Duality  pairings  between  spaces  and  their 
duab  will  be  denoted  by  (•,  ■}.  Norms  of  functions  belonging  to  and  H'( r)  are 

denoted  by  ||  •  ||,  and  |)  •  ij^p,  respectively.  Of  particular  interest  are  the  L2(Q)-norm 
||  •  ||o,  the  norm 

iwb-eI^-I’+iwb. 

j=t  II  oxi  llo 

and  the  f/2(fi)-norm 


Corresponding  Sobolev  spaces  of  real,  vector- valued  functions  having  r  compo¬ 
nents  will  be  denoted  by  H*(2>),  e.g.,  H^D)  =  [Hl(V)]r.  Of  particular  interest  will 
be  the  spaces  L2(D)  =  H°(Z>)  =  [L3(D)]r, 

H\V)  =  {Bj  G  L\V)  |  |SL  e  L\V)  for  j  =  1 . r  and  k  =  1 . dj  , 


and 


" 'm = {•» « I  *  L'm  ■ 

for  j  =  1 , . . . ,  r  and  k,l=  1 , . . . ,  d  j  , 

where  vj,  j  =  1, . . . ,  r,  denote  the  components  of  v.  Also  of  interest  is  the  subspace 

Hj(D)  =  {v€H1(P)  |  =  0  on  dV,  j  =  1, . . . ,  r  |  . 

Norms  for  spaces  of  vector-valued  functions  will  be  denoted  by  the  same  notation  as 
that  used  for  their  scalar  counterparts.  For  example, 

IMI?  =  £lM?  and  IMI*  r  =  ]£  INI*  r  • 

j= i  ;= i 

We  denote  the  L2(fi)  and  L2(fl)  inner  products  by  (•,  •),  i.e.,  for  p,  q  G  L2(Q)  and 
u,  v  G  L2(D) 

(Pi  ?)  =  I  pqdQ  and  (u,  v)  =  /  u  •  v  dQ . 

Jn  Jn 

Similarly,  we  denote  by  (•,  -)p  the  L2(T)  and  L2(r)  inner  products,  i.e.,  for  p,  q  G  L2(T) 
and  u,v  G  La(r) 


and 


u  •  v  dr . 
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Since  in  all  cases  L3-spaces  will  be  used  as  pivot  spaces,  the  above  inner  product 
notation  can  also  be  used  to  denote  duality  pairings  between  functions  defined  on 
H‘- spaces  and  their  dual  spaces. 

For  details  concerning  the  notation  employed,  one  may  consult,  e.g.,  [1]. 

4.1.  Distributed  controls  for  the  von  Kdrmdn  plate  equations 


For  this  application  we  will  use  distributed  controls,  i.e.,  control  is  effected  through 
a  source  term  in  the  governing  partial  differential  equations.  Let  0  be  a  bounded,  con¬ 
vex  polygonal  domain  in  R2  and  let  T  denote  the  boundary  of  ft.  The  von  Karman 
equations  for  a  clamped  plate  are  given  by  (see,  e.g.,  [9]  or  [18]) 


A2^i  +  =  0 


in  ft , 


A2t/>2  -  [tp\ .  f/’a]  =  A g  in  ft , 


and 


where 


dtp  i  dtpi 


on  r , 


r/  „  .  P± d24>  ,  dH d*<p  n  d2tp  d'<p 
'  dx\  dx\  dx\  $x%  dx\x-i  dx\x2 


Here,  tp\  denotes  the  Airy  stress  function,  tp2  the  deflection  of  the  plate  in  the  direction 
normal  to  the  plate,  A g  is  an  external  load  normal  to  the  plate  which  depends  on  the 
loading  parameter  A,  and  d(-)/dn  the  normal  derivative  in  the  direction  of  the  outer 
normal  to  I\ 

By  introducing  appropriate  rescalings,  i.e.,  by  replacing  tpi  by  A^i,  tp2  by  Xip2, 
and  g  by  Ay,  we  can  rewrite  the  von  Karman  equations  as  follows: 


(4.1) 


A2V>1  +  •ztylytpl}  =  0 


in  ft , 


(4.2) 
and 

(4.3) 


A2V>2  —  A  [0i ,  ^2]  =  A  g  in  ft 


dtpi  dtp2 


on  r. 


We  introduce  the  spaces 

ffo2(«)  =  {v> e  tf2(ft)  |  tp  =  0,  ^  =  0  on  rj  , 

Hg(ft)  =  [ff2(ft)]2,  H~ *(ft)  =  (//02(ft))*  ,  and  H"2(ft)  =  (Hg(ft))* 
and  the  bilinear  form 


a(tP,<P)  =  /  AxpA<pdtl  V  tp,<p  €  //2(ft) 

Ja 
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in  order  to  define  the  following  weak  formulation  of  the  von  Karman  equations  (4.1)- 

(4.3) :  find  4  =  (tj>i ,  V'a)  €  Hjj(O)  such  that 

(4.4)  «(*,*i)+  ^(totojto)  =  0  V^i  €  //02(fl) 

and 

(4.5)  atoto)  -  A((V»i,^a],^a)  =  k{g,<h)  V  ds  €  ^o(^)  • 

Using  the  identity 

(4.6)  (to  4],  C)  =  (14,  a  4)  V  4, 4,  <  e  //02(n) , 

one  can  show  that  for  each  g  €  H~7(Q),  (4.4)-(4.5)  possesses  at  least  one  solution 
4  =  (4i  1 4i)  €  Hj(fi)  and  that  all  solutions  of  (4.4)-(4.5)  satisfy  the  a  priori  estimate 

(4-7)  |to||3  +  M2<C||ff||_2; 

see,  e.g.,  [18],  for  details.  In  the  sequel  a  solution  to  (4.1)-(4.3)  will  be  understood  in 
the  sense  of  (4.4)-(4.5). 

Given  a  desired  state  4o  —  (V’lOi 4v>)  €  L2(Q),  we  define  for  any  4  —  (4\,4l)  € 
H3(0)  and  g  G  LJ(Q)  the  functional 

J(4,g)  =  J(4\,4%,g) 

(4‘8)  =  \  f  (to  -  4io)7  +  to  -  Iho)2)  dft  +  \  jf  g7  dQ . 

We  then  consider  the  following  optimal  control  problem  associated  with  the  von 
Karman  plate  equations: 

(4.9)  min  {  J(4,  g)  I  4  €  H2(fl) ,  g  G  0  }  subject  to  (4.4)-(4.ti) , 
where  6  is  a  subset  of  L3( fl). 

We  define  the  spaces  X  =  Hg(O),  Y  =  H -3(fl),  G  =  L3(fl),  and  Z  =  V{(1). 
By  compact  imbedding  results,  Z  «-+«-♦  Y.  For  the  time  being,  we  assume  that  the 
admissible  set  6  for  the  control  g  is  a  closed,  convex  subset  of  G  =  L7(U). 

Let  the  continuous  linear  operator  T  G  C(Y ;  X)  be  defined  as  follows:  for  f  = 
(fi  ,fi)£Y  —  H~3(fi),  Tf  =  4  €  X  =  H3(fi)  is  the  unique  solution  of 

<*toto)  =  (fi,4i)  V  4i  €  H7(Sl) 

and 

<*to>  4i)  =  (fi,4 2)  V  fa  G  H*(Q) . 

It  can  be  easily  verified  that  T  is  self-adjoint. 

We  define  the  (differentiable)  nonlinear  mapping  N:  X  —>Y  by 
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or  equivalently 


(W(0),0)  =  j  ([0a.0a].0i)  -  ([0i,0a]>0a)  V  0  =  (0i>0a)  €  X 
and  define  K  :  g  €  £3(fl)  — ►  Y  by 

*~C)- 

or  equivalently, 

{Kg  A)  =  -(g,h)  V*  =  (*i,*a)€X. 

Clearly,  the  constraint  equations  (4.4)-(4.5)  can  be  expressed  as 

0  +  A7W(0)  +  ATKg  =  O, 

i.e.,  in  the  form  (2.2).  With  the  obvious  definitions  for  )  and  £{■),  i.e., 

r(4)  =  \Jn ({1>x  -  tho)2  +  (02 - 0ao)3) dti  v 0 ex 

and 

£{a)  =  ^ Jn92 dil  Vs e G, 
the  functional  (4.8)  can  be  expressed  as 

«7(0,*)  =  A:F(0)  +  A  £{g), 

i.e.,  in  the  form  (2.1).  Thus,  the  minimization  problem  (4.9)  is  in  the  form  of  the 
minimization  problem  (2.3). 

We  are  now  in  a  position  to  verify,  for  the  minimization  problem  (4.9),  all  the 
hypotheses  of  §2  and  §3. 

4.1.1.  Verification  of  the  hypotheses  for  the  existence  of  optimal  solutions. 

We  first  verify  that  the  hypotheses  (H1)-(H6)  hold  in  the  current  setting. 

(HI)  is  obviously  satisfied  with  a  lower  bound  0. 

(H2)  holds  with  a  =  1  and  /?  =  2. 

(H3)  is  verified  with  the  choice  (0^,0^)  €  X  x  0,  where  gW  is  an  arbitrarily 
chosen  element  in  0  and  0^  =  (0^ ,  t/4°^)  is  a  solution  of 

^ V(10)  +  |[02O) .  02O)]  =  0  in  n . 


A20<o)-A[0<o>,4o>]  =  A?<°>  infl, 


In  order  to  verify  (H4),  we  assume  {$[(")}  C  0  is  a  sequence  satisfying  — *  g 

in  L2( n);  then,  we  have  g^  — -  g  in  so  that  lim„_oo(p^n\  z)  —  {g,z)  for 

all  z  €  W3(n),  i.e.,  Kg^  —  Kg  in  Y.  Assume  that  the  sequence  {0^}  C  H§(Q) 


satisfies  *(n)  -  0  in  then,  {d2+(n) / tondtj)  -  {d2+/dzidxj)  in  Ls(ft)  and, 

by  using  a  compact  imbedding  result,  4^  — »  0  in  L2(0).  Now,  using  the  identity 
(4.6), 

(N(^),4)  =  ([0<n),0<n)],0,)  -  ([V(1B),^n)],^) 

=  |  ([V4n),<H,t4n))  -  ([V'in),^2]»  V'a”>) 

-*  ^  {[03,  0i],  0a)  -  (E01 . 03],  V»a) 

=  ^  ([03, 02],  0l)  “  ([01, 03],  03)  =  (#(0),#)  • 

Hence,  (H4)  is  verified. 

The  verification  of  (H5)  follows  directly  from  the  observation  that  the  mappings 
4  ■_  T(4)  =  (1/2)H0  -  0o||B  and  ~  £(g)  =  (l/2)||ff||g  are  convex. 

The  verification  of  (H6)  is  a  trivial  consequence  of  the  a  priori  estimate  (4.7). 

It  is  now  just  a  matter  of  citing  Theorem  2.1  to  prove  the  existence  of  an  optimal 
solution  that  minimizes  (4.8)  subject  to  (4.4)-(4.5). 

THEOREM  4.1.  There  exists  a  (4,9)  €  H2(fi)  x  0  such  that  (4.8)  is  minimized  subject 
to  (4.4)-(4.5).  0 

4.1.2.  Verification  of  the  hypotheses  for  the  existence  of  Lagrange  multipliers. 

We  now  assume  (4,  g)  is  an  optimal  solution  and  turn  to  the  verification  of  hy¬ 
potheses  (H7)-(H9). 

The  validity  of  (H7)  is  obvious. 

(H8)  holds  since  the  mapping  g  >-*  £(g)  =  (l/2)||<y||o  is  convex. 

(H9)  can  be  verified  as  follows.  For  any  4  G  X,  the  operator  N'(4)  ■  X  — <•  Y  is 
given  by 

Thus,  using  the  definition  of  [•,  •],  we  obtain  that  N'(4)  •  4  G  LJ(ft)  =  Z. 

The  Lagrangian  is  given  by 

m>,9,r),k)  =  *  {a{4i,m)+  ^([02, 0s],  »7i) 

+  a(42,  m)  -  ^([0i,  02],  m)  -  Hg,  »ft)} 

for  all  (0,  Ar)  eXxGxXxR  =  Hg(fi)  x  L2(fi)  x  Hg(fi)  x  R.  Note  that  in  this 
form  of  the  Lagrangian,  the  Lagrange  multiplier  t/  €  X  =  Y*  so  that  we  have  already 
introduced  the  change  of  variables  indicated  between  (2. 17)-(2. 18)  and  (2.19)-(2.21). 

Having  verified  the  hypotheses  (H7)-(H9),  we  may  apply  Theorem  2.4  to  conclude 
that  there  exists  a  Lagrange  multiplier  ij  £  X  =  H2(Q)  and  a  real  number  k  such  that 

(4.10)  11  +  XT* ([JV'(0)]*  r,~  kJ^(4,9))  =  0 
and 

(4.11)  -C(0,  flr,  *)  <  £(0, 3r,  *)  V  z  6  0 


and  that  for  almost  all  values  of  A,  we  may  choose  k  =  1. 


Recall  that  T  is  self-adjoint.  Also,  note  that  for  any  G  X  =  H2(ft), 

Thus,  (4.10),  with  k  =  1,  can  be  rewritten  as 

(4.12)  a(Ci.ili)  -A([^,  vjJ.Ci)  =  A(V>,  -  fco.Ci)  V(,  G  f/02(fl) 

and 

14  131  a(C2.»U)  +  A((V'3,»;i],C3)  -  AQth.iftJ.Cs) 

'■  '  =A(^-V>30,C3)  VC2  €f/02(ft). 

Using  the  definition  of  the  Lagrangian  functional,  (4.11),  with  k  =  1,  can  be  rewritten 
as 

|  (z,  z)  +  A(z,  I?a)  -  ^  (y.  ff)  -  A(ff ,  r&)  >  0  V  z  G  0 . 

Note  that,  in  the  above  expression,  we  have  already  employed  hypothesis  (H12)  which 
in  the  current  context  is  trivially  satisfied  with  E  the  identity  operator  on  G*  =  G  — 
L2( Q).  For  each  t  €  (0, 1)  and  each  t  G  0,  set  z  =  rt  +  (l -t)g  G  0  in  the  last  equation 
to  obtain 

eJ 

y  (*  ~  9,  *  -  ff)  +  *  (1  ~  9, 9)  +  ( (t  -  g,  in)  >  0  V  t  G  0 
so  that,  after  dividing  by  i  >  0  and  then  letting  c  — ►  0+ ,  we  obtain 

(4.14)  (t~g,g  +  m)>0  vie©. 

We  see  that  for  almost  all  values  of  A,  necessary  conditions  for  an  optimum  are 
that  (4.4)-(4.5)  and  (4.12)-(4.14)  are  satisfied.  The  system  formed  by  these  equations 
will  be  called  an  optimality  system. 

We  now  specialize  to  the  case  0  =  L2(Q).  Note  that  the  hypothesis  (H 10)  is 
satisfied.  Then,  using  Theorem  2.5,  we  see  that  the  inequality  (4.14)  becomes  an 
equality  and,  by  letting  z  —  t  —  g  vary  arbitrarily  in  L2(Sl),  we  now  have,  instead  of 

(4.14) , 

(4.15)  (z,y  +  i|  3)  =  0  VzGL2(«). 

Thus,  according  to  that  theorem,  we  have  that  for  almost  all  A,  an  optimality  system 
of  equations  is  now  given  by  (4.4)-(4.5),  (4.12)-(4.13),  and  (4.15).  However,  we  can  go 
further  and  verify  that  the  hypothesis  (Hll)'  is  valid,  which  in  turn  will  justify  the 
existence  of  a  Lagrange  multiplier  satisfying  the  optimality  system  for  all  A  G  A.  We 
now  assume  the  domain  Q  is  a  convex  polygon  with  no  angles  greater  than  126°. 

Let  A  be  given  such  that  1  /A  is  an  eigenvalue  of  —TN'(i>),  where  (f/>,  g)  6  H2(il)  x 
L2(fi)  is  an  optimal  pair  that  minimizes  (4.8)  subject  to  (4.4)-(4.5).  We  wish  to  show 
that  for  each  f  €  H-2(J1),  there  exists  a  g  G  L2(fi)  and  a  ^  G  H2(12)  such  that 

*  +  A7W'(*)-*+A:TRy  =  f, 


i.e., 

(4.16)  o(^,0,)  +  A([V>2,^2]^i)  =  (/i^i)  V  G  Hq(Q.) 
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and 

-  -Ms.**) 


V^€W0J(«)- 

The  existence  of  such  a  *  can  be  shown  in  a  manner  similar  to  that  for  showing  the 
existence  of  a  solution  to  the  von  Karman  equation;  the  key  step  is  that  by  adding  the 
two  equations  with  the  test  function  *  replaced  by  *,  we  have  the  a  priori  estimate 

o(*i ,  V>i)  +  a{i>3,i>i)  =  {fi,4i)  +  (/a>  V>a)  ■ 

Then,  we  choose  g  =  —  [Vq,  *2].  Note  that  regularity  results  for  the  biharmonic  equar- 
tion  applied  to  (4.4)-(4.5)  yield  4  €  H4(fi)  (see  [3]).  Hence,  using  imbedding  theorems 
we  deduce  that  g  £  L2(fl).  It  is  obvious  that  g  and  4  satisfy  (4.16)-(4.17),  i.e.,  we  have 
verified  (Hll)'.  Hence  we  conclude  that  for  all  A,  the  optimality  system  (4.4)-(4.5), 
(4.12)-(4.13),  and  (4.15)  has  a  solution.  Thus,  we  have  Theorem  2.6  which,  in  the 
present  context,  is  given  as  follows. 

THEOREM  4.2.  Lei  (4,  g)  £  H2(fl)  x  L3(Q)  denote  an  optimal  solution  that  minimizes 
(4.8)  subject  to  (4.5)-(4.6).  Then,  for  all  A  6  A,  there  exists  a  nonzero  Lagrange 
multiplier  i|  £  H2(Q)  satisfying  the  Euler  equations  (4.12)-(4.13)  and  (4.15).  □ 

4-1.3.  Verification  of  the  hypotheses  for  approximations  and  error  estimates. 

We  finally  verify  the  hypotheses  (H13)-(H19)  that  are  used  in  connection  with 
approximations  and  error  estimates. 

A  finite  element  discretization  of  the  optimality  system  (4.4)-(4.5),  (4.12)-(4.13), 
and  (4.15)  is  defined  in  the  usual  manner.  We  first  choose  families  of  finite  dimensional 
subspaces  Xh  C  Hj|(Q)  and  Gh  C  L2(fi)  parameterized  by  a  parameter  h  that  tends 
to  zero  and  satisfying  the  following  approximation  properties:  there  exists  a  constant 
C  and  an  integer  r  such  that 

(4.18)  inf  !|*  -  *fc||2  <  CW"|Wlm+a  -  V  *  G  Hm+2(fi) ,  1  <  m  <  r 
and 

(4.19)  inf  ||*  -  zh||0  <  Chm||z||m  ,  V  z  £  Hm(Q) ,  1  <  m  <  r . 

xk£Gk 

One  may  consult,  e.g.,  [8]  for  some  finite  element  spaces  satisfying  (4.18)  and  (4.19). 
For  example,  one  may  choose  Xh  =  V*  x  V*  where  Vh  is  the  piecewise  quintic-C'(Q) 
finite  element  space  constrained  to  satisfy  the  given  boundary  conditions  and  defined 
with  respect  to  a  family  of  triangulations  of  Q.  In  this  case,  h  is  a  measure  of  the  grid 
size.  For  simplicity,  one  may  choose  Gh  =  Vh. 


(4.17) 


a(*a,  4i)  ~  A([*i .  *a],  4t)  ~  A([*t  >  &]>  ^2) 


To  show  this,  we  first  let  *  6  H2(0)  be  a  solution  of 

<»(*i.*i)  +  =  (/tt*i) 

and 

a(*j.*a)  -  —  (h>4  a) 
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Once  the  approximating  spaces  have  been  chosen,  we  may  formulate  the  approx¬ 
imate  problem  for  the  optimality  system  (4.4)-(4.5),  (4. 12)-(4. 13),  and  (4.15):  seek 
+k  £  Xk ,  gk  £  Gk ,  and  tf  £  Xk  such  that 

(4.20)  ««**)+ 0  V  £Vk, 

(4.21)  a(V»J,^)  -  X([rJ>k,  tf],  <f>k)  =  (g,  4>k)  V  4>k  £  l/h  , 

(4.22)  a(C,k,i?f)  -  A([V>§,  =  A(V<?  -  ^-io,Ci*)  V  Cf  €  Vk  , 

(4.23)  a(Cah,»J2k)  +  A([^,^],C2fc)  -A([tf, «£].<$)  =  A(</>k  -  02O,C2h)  Vtf  €  V*  , 

and 

(4.24)  (*\  #*  +  •£)«  0  V:k€Gk. 

The  operator  Th  €  £(Y ;  Xk)  is  defined  as  follows:  for  f  £  Y ,  Th f  =  tf>h  £  Xh  is 
the  unique  solution  of 


«(*?,*? )  =  (/t,tf> 

and 

a(V-$,^)  =  (/a,^)  V^€Kk. 

Since  T  =  71*,  we  define  (T*)k  =  T*. 

We  define  the  operator  £k  :  L2(ft)  — ♦  Gh  as  the  L2(ft)-projection  on  Gk,  i.e.,  for 
each  j  £  L2(ft), 

(Ehg,^)  =  (g,<t>h)  V^6Ck. 

Since  G  =  I2  (ft)  is  reflexive,  Eh  is  in  fact  an  operator  from  G *  Gh. 

By  the  well-known  results  concerning  the  approximation  of  biharmonic  equations 
(see,  e.g.,  [2]  or  [8]),  we  obtain 


||(T-Tk)f||x-0 

as  h  -*  0,  for  all  f  €  Y.  This  is  simply  a  restatement  of  (H13). 

(H14)  follows  trivially  from  (HI 3)  and  the  fact  that  T  is  self-adjoint  and  we  have 
chosen  (T*)k  =  Th. 

(H15)  follows  from  the  best  approximation  property  of  L2(ft)-projections  and 
(4.19). 

(H16)  and  (H17)  follow  from  the  fact  that  N  and  T  are  polynomials.  Here  we 
also  use  imbedding  theorems  and  Cauchy  inequalities. 

We  set  Z  =  Z  =  L*(ft).  For  each  ij  6  H2(ft)  and  C  €  H2(ft),  Sobolev  imbedding 
theorems  imply  that 


M*-?= 


(  -[V’s.’fe]  \ 

V[V>2.»;i]  -  bPum]J 


€Z, 


23 


and 

(rwt),=  te)d' 

These  relations  verify  (H18). 

From  the  deiinition  of  the  operator  A1  we  see  that  K  maps  L2(fi)  into  L](n),  i.e., 
K  maps  G  into  Z.  Thus  (H19)  is  verified. 

Hence,  we  are  now  in  a  position  to  apply  Theorem  3.5  to  derive  error  estimates 
for  the  approximate  solutions  of  the  optimality  system  (4.4)-(4.5),  (4. 1 2)-(4. 13)  and 
(4.15).  It  should  be  noted  that  Lemma  3.4  implies  that  for  almost  ail  values  of  A,  the 
solutions  of  the  optimality  system  are  regular. 

THEOREM  4.3.  Assume  that  A  is  a  compact  interval  of  R+  and  that  there  exists 
a  branch  {(A,Vl(A),j/(A),i|(A))  :  A  €  A }  of  regular  solutions  of  the  optimality  system 
(4.4)-(4.5),  (4.12)-(4.13),  and  (4.15).  Assume  that  the  finite  element  spaces  Xh  and 
Gh  satisfy  the  hypotheses  (4. 18)-(4. 19).  Then,  there  exists  a  6  >  0  and  an  h0  >  0  such 
that  for  h  <  ho,  the  discrete  optimality  system  (4.20)-(4.24)  has  a  unique  branch  of 
solutions  {(A,lfrk(A),<jfe(A),9*(A))  :  A  6  A)  satisfying 

{||*\A)  -  *(A)lh  +  \\9hW  ~  ff(A)||o  +  IN* (A)  -  «J(A)||2)  <  6  for  all  A  G  A. 
Moreover, 

lim{||*k(A)  -  *(A)||2  +  ||/(A)  -  s(A)||0  +  \\r,h(\)  -  9(A)||2}  =  0 , 
uniformly  in  A  G  A. 

If,  in  addition,  the  solution  of  the  optimality  system  satisfies  (V’(A),  g(A),  »|( A))  G 
Hm+2(ft)  x  Hm({2)  x  Hm+2(il)  forX  G  A,  then  there  exists  a  constant  C,  independent 
ofh,  such  that 

IMA)  -  *h(A)||2  +  ||s(A)  -  /(A)||o  +  ||9(A)  -  9h(A)||2 
<  C/.m(|MA)||m+2  +  ||,(A)||m  +  ||9(A)||m+2)  , 

uniformly  in  A  G  A. 

Proof :  All  results  follow  from  Theorem  3.5.  For  the  last  result,  we  also  use  (3.25)  and 
the  estimates  (see,  e.g.,  [2]  or  [8]) 

||(ThT-1  -  I)1>fo  <  C/.m|M|m+2  for  *  G  Hm+2(fl) , 

||((T*)h(T*)-1  -  /)9||2  =  ||(TkT-1  -  /)i,||2  <  Chm\\n\\m+2  for  9  G  H”>+2(Q) , 

and 

-  /)j||0  <  Chm\\g\\m  for  g  G  . 

In  these  estimates,  the  constant  C  is  independent  of  h,  if>,  g,  9,  and  A.  □ 

Remark.  In  fact,  we  obtain  from  (4.15)  that  g  =  —  r?2  so  that  the  term  ||p(A)j|m  in  the 
right-hand  side  of  the  error  estimate  is  redundant.  0 

Remark.  By  using  (4.15)  again, .along  with  (4.24)  and  the  error  estimate  in  Theorem 
4.3,  we  have  the  following  improved  error  estimate  for  the  approximation  of  the  control 

ll«(A)  -  /(A)||2  =  1193(A)  -  9$(A)||a  <  Chm(|MA)||m+2  +  l«A)||m+a)  • 
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Of  course,  we  also  use  the  fact  that  we  have  chosen  Gk  —  Vk  C  H2(Sl).  □ 

4.2.  Neumann  boundary  controls  for  the  Ginzburg* Landau  superconduc¬ 
tivity  equations 

For  this  application  we  will  use  Neumann  boundary  controls,  i.e.,  control  is  ef¬ 
fected  through  the  data  in  a  Neumann  boundary  condition.  Let  0  be  a  bounded  open 
domain  in  R*,  d  —  2  or  3,  and  let  T  be  its  boundary.  A  simplified  Ginzburg-Landau 
model  for  superconductivity  is  given  by 


-A^i  +  (V> 2  +  tp 2  +  |A|2  -  l)  tpi  -  V  •  (A fa)  -  A  •  Vt/>2  =  0  in  0 , 


and 


-Ai>2  +  (V>1  -¥'P 2  +  |A|2  -  l)  rp2  +  v  •  (AV»i)  +  A  ■  V^i  =0  in  Q , 
n  •  (VV<i  +  AV>2)  =  on  T , 

n  •  (V^>2  -  A^i)  =  Xg2  on  f. 


Here,  Vi  and  -02  denote  the  real  and  imaginary  parts,  respectively,  of  the  complex¬ 
valued  order  parameter,  A  is  a  given  real  magnetic  potential,  g i  and  g2  are  related 
to  the  normal  component  of  the  current  at  the  boundary,  and  A  >  0  is  a  “current 
loading”  parameter.  These  equations  are  a  special  case  of  a  more  general  model  for 
superconductivity  wherein  A  is  also  unknown;  see,  e.g.,  [22]  for  a  derivation  of  the 
general  model.  It  can  be  shown  that  in  certain  limits,  e.g,  high  values  of  the  applied 
field,  the  above  simpler  model  is  valid;  see  [7]. 

By  introducing  appropriate  rescalings,  i.e,  by  replacing  rpj  by  A rpj  and  gj  by  A g,, 
j  =  1,2,  we  can  rewrite  the  above  Ginzburg-Landau  equations  as  follows: 


(4.25) 

-AV>i  +  (|A|2- 

(A^>2)  -  A  •  Vtp 2  +  A 

{'Pi  +  'Pi)  ^1=0 

in  0 , 

(4.26) 

-AV>2  +  (|A|2  - 

1 

i- 

+ 

<3 

(A^i)  +  A  •  -(-  A 

{'Pi  +  'Pi)  ^2  =  0 

in  SI , 

(4.27) 

u  -  (VVq 

+  Atp2)  =  Xg\  on  T 

) 

and 

(4.28) 

n  •  (V02 

-  Atp\)  =  A <72  on  f 

We  introduce  the  bilinear  forms 


*(*, <f>)  =  J(v+  ■  V*  +  (|A|2  -  I)**)  dfi  V  €  H'(Sl) 


and 


b(r/>,</>)=  f  A  ■  (rpv<p  -  <pVrp)  dU  V  rp,  <P  €  Hl{Sl) . 
J  n 

We  assume  that  A  €  H*(fi).  Note  that 

a (rj>,<p)  =  a(<p,rp)  and  6(0, 0)  =  -b(<p,rp) . 
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Then,  a  weak  formulation  of  the  Ginzburg-Landau  equations  (4.25)-(4.28)  is  defined 
as  follows:  seek  $  =  (V»i ,  ^2)  €  H1  (Q)  such  that 

(4.29)  a{ii>i,<l>i)  +  b{fa,<j>\)  +  =  H9i,<t>i)r  V  4>\  £Hl{&) 


(4.30)  a(fa,fa)  -  b(fa,fa)  +  A((V>]  +  *l>\)fa,fa)  =  ^(92,fa)r  V  fa  £  H J(f2) . 

It  can  be  shown  that,  for  each  g  =  ($1,02)  £  H-1/2(r),  (4.29)  and  (4.30)  possess  at 
least  one  solution  if>  £  H*(Q)  and  that  all  solutions  of  (4.29)  and  (4.30)  satisfy  the  a 
priori  estimate 


(4.31) 


IIV’illi  +  IhMIi  <  C  (llst||-i/2,r  +  Ilff2||-i/2,r)  ; 


see,  e.g.,  [11],  for  details.  In  the  sequel,  a  solution  of  (4.25)-(4.28)  will  be  understood 
in  the  sense  of  (4.29)-(4.30). 

Given  a  desired  state  i>0  —  (fao,  ^20)  £  L2(fi),  we  define  for  any  tj>  =  (\pi ,  fa)  £ 
H^fi)  and  g  =  (31,02)  G  L2(r)  the  functional 

(4.32)  J{fa  g)  =  \  J  ({fa  -  V>io)2  +  {fa  -  fao)2)  dfi +  £  J (,?  +  022)  dr . 

We  then  consider  the  following  optimal  control  problem  associated  with  the  Ginzburg- 
Landau  equations  for  superconductivity: 

(4.33)  min  (  J(fa  g)  |  1>  6  H^fl) ,  g  6  0  }  subject  to  (4.29)  and  (4.30) , 
where  0  is  a  subset  of  L2(r). 

We  define  the  spaces  X  =  H^ft),  Y  =  (H1^))*,  G  =  L2(r),  and  Z  - 
[Hi/2+‘(n)]«  where  e  €  (0, 1/2)  is  chosen  such  that  Hl(G)  H1/2+e(G) «-»  L4(ft). 
By  compact  imbedding  results,  L4/,3(f2) «-»  Z  Y.  For  the  time  being,  we  assume 
that  the  admissible  set  0  for  the  control  g  is  a  closed  convex  subset  of  G  =  L2(T). 

Let  the  continuous  linear  operator  T  £  £(Y;  X)  be  defined  as  follows:  for  each 
f  =  {fufi)  €  V  =  (H1  («))*,  Tf  =  £  X  =  H^Il)  is  the  unique  solution  of 

<*{*l>i  >  fa)  +  b{fa,  <f>i)  =  V  <pi  £  H^fi) 


a{fa,  ^2)  -  b{fa ,  fa)  =  {/2,  fa)  V  fa  £  H1( fi) . 

It  can  be  easily  verified  that  T  is  self-adjoint.  Also,  it  can  be  shown  that  for  most 
choices  of  A,  the  operator  T  is  well  defined;  see  [11]. 

We  define  the  (differentiable)  nonlinear  mapping  X  :  X  — *  Y  by 


N(j,)  _  (Wi  +  +Dfa  ^ 


v^a 


or  equivalently 


(JV {fa, fa  =  ((V*2  +  ^2)^1,  <f>i)  +  ((^2  +  fal)fa,fa)  V  4  =  {fa,  fa)  €  X 
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A;J\  ■  ■ 


and  define  K  :  H-‘/J(r)  -  Y  as  the  injection  mapping: 

(Ku,  v)  =  -(a,  v)r  V  t  €  H'^r)  ,V  v  6  H^Q) . 

Clearly,  the  constraint  equations  (4.29)-(4.30)  can  be  expressed  as 

+  +  \TN(+)  +  XTKg  =  0, 

i.e.,  in  the  form  (2.2).  With  the  obvious  definitions  for  T(-)  and  £(•),  i.e., 

r(4)  =  \  jf  ((* x  -  tMJ  +  (*a  -  *1  o)2)  dQ  v * e  x 

and 

€{g)=\jyx+9l)dT  Vg  GG, 
the  functional  (4.32)  can  be  expressed  as 

J(*,9)  =  A:F(*)  +  A£(g), 

i.e.,  in  the  form  (2.1).  Thus,  the  minimization  problem  (4.33)  is  in  the  form  of  the 
minimization  problem  (2.3). 

We  are  now  in  a  position  to  verify,  for  the  minimization  problem  (4.33),  all  the 
hypotheses  of  §2  and  §3. 

4-2.1.  Verification  of  the  hypotheses  for  the  existence  of  optimal  solutions. 

We  first  verify  that  the  hypotheses  (HI)-(H6)  hold  in  the  current  setting. 

(HI)  is  obviously  satisfied  with  a  lower  bound  0. 

(H2)  holds  with  a  =  1  and  fi  =  2. 

(H3)  is  verified  since  *  =  0  and  g  =  0  is  obviously  a  solution  of  (4.29)-(4.30). 

In  order  to  verify  (H4),  we  assume  {g^}  C  0  C  L2(."'  is  a  sequence  satisfying 
g(n)  —  g  in  L2(r);  then,  we  have  g(n)  — *  g  in  H-1^2(r)  so  nat  limn_0„(g*n\  v)r  = 
(g.v)r  for  all  v  e  H1(f2),  i.e.,  KgW  K g  in  Y .  Assume  that  the  sequence 

{*^n  }  C  H*(fi)  satisfies  — -  *  in  H1(fl);  then,  by  using  the  compact  imbed¬ 

ding  H‘(0)  «-♦«-»  L4(fi),  *(n)  -  *  in  L4(Q).  Now, 

(n^a) = (((v-(r))2 + +  (((4n))2 + 

-*•  ((*?  +  *|)lh>*i)  +  ((*i  +^2)^2,^)  =  {N{1>),4>). 

Hence,  (H4)  is  verified. 

The  verification  of  (H5)  follows  directly  from  the  observation  that  the  mappings 
4  >-*  T{4)  =  (1/2)||*  -  *ollo  and  g  *-*■  £(g)  =  (l/2)||g|lo,r  are  convex. 

The  verification  of  (H6)  is  a  trivial  consequence  of  the  a  priori  estimate  (4.31). 

It  is  now  just  a  matter  of  citing  Theorem  2.1  to  prove  the  existence  of  an  optimal 
solution  that  minimizes  (4.32)  subject  to  (4.29)-(4.30). 

THEOREM  4.4.  There  exists  a  (*,  g)  €  such  that  (4.32)  is  minimized  subject 

to  (4.29)-(4.30).  □ 

4-2.2.  Verification  of  the  hypotheses  for  the  existence  of  Lagrange  multipliers. 
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We  now  assume  (^,g)  is  an  optimal  solution  and  turn  to  the  verification  of 
hypotheses  (H7)-(H9). 

The  validity  of  (H7)  is  obvious. 

(H8)  holds  since  the  mapping  g  £ (g)  =  (1/2)  fr  |g|2  dr  is  convex. 

(H9)  can  be  verified  as  follows.  For  any  ^  €  X,  the  operator  i V'(^)  :  A'  — ►  V'  is 
given  by 

A  —  (  ^2)^1  +  (2^1  ^2)^2  \  W  A  —  (A  A  \  e  X 

N  (,t) '  *  -  (  (3#|  +  +  (2*,  fe  )*, ) 

Thus,  we  obtain  that  TV'(^)  •  +  6  L“'3(fi)  --  [H1/2+‘(fi)]*  =  Z. 

The  Lagrangian  is  given  by 

C(+,g,tl,k)  =4J(*,  g) 

-  {a(^i.»?i)  +  +  ^2)^1.  Vi)  ~  A(yi,»?i)r 

+  a(^2,»j2)  -  ,<j>2)  -  A((V>?  +  T)i)  -  A(g2,»)2)r} 

for  all  (^,  g,  Jfc)  6  X  x  G  x  X  x  R  =  H1  (Q)  x  L2(r)  x  H1  (ft)  x  R .  Note  that  in  this 
form  of  the  Lagrangian,  the  Lagrange  multiplier  ij  6  X  =  Y’  so  that  we  have  already 
introduced  the  change  of  variables  indicated  between  (2. 17)-(2. 18)  and  (2.19)-(2.21). 

Having  verified  the  hypotheses  (H7)-(H9),  we  may  apply  Theorem  2.4  to  conclude 
that  there  exists  a  Lagrange  multiplier  q  €  X  =  H^ft)  and  a  real  number  k  such  that 


(4.34) 
and 

(4.35) 


*1  +  A  r  ([W W  t,-kj^,  g))  =  0 


£(+,g,*i,k)<c(+,z,ii,k)  v*ee 


and  that  for  almost  all  values  of  A,  we  may  choose  4=1. 

Recall  that  T  is  self-adjoint.  Also,  note  that  for  any  i/>  €  X  —  H](ft), 

[N  <*>J  •  ■ "  -  ( (3Vi  +  *?)-«  +  (»|  Ml, )  v  1  -(->'.  *)  £  X. 

Thus,  is  self-adjoint  as  well  and  (4.34),  with  4=1,  can  be  rewritten  as 

«(Ci .  m)  ~  KCi >  %)  +  A((3t/>i  +  Ci) 


(4.36) 
and 

(4.37) 


+  A((2^i^2)»?2>Ci)  =  A(V>i  -  ^10,6)  V  Ci  Gf/'(ft) 

a(C2.»?2)  +  +  A((3V|  +  )^2,  C2) 

+  A  ((2V»1  >  C2)  =  A(V>2  -  V,,2o>  C2)  V  C2  €  //'(ft)  ■ 


Using  the  definition  of  the  Lagrangian  functional,  (4.35),  with  4=1,  can  be  rewritten 
as 

^(z,z)r  4-  A(z,ij)r  -  ^(g,g)r-  A(g,ij)r  >0  V  z  e  @  • 
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Note  that,  in  the  above  expression,  we  have  already  employed  hypothesis  (B12)  which 
in  the  current  context  is  trivially  satisfied  with  E  the  identity  operator  on  G*  =  G  = 
LJ(r).  For  each  f  6  (0, 1)  and  each  t  €  0,  set  a  =  ct  +  (1  -  f)g  €  0  in  the  last 
equation  to  obtain 

y  (*  -  g.  t  -  g)r  +  <  (t  -  g,  g)r  +  f  (t  -  g,  »?)r  >  0  V  t  G  © 

so  that,  after  dividing  by  t  >  0  and  then  letting  f  — *  0+,  we  obtain 

(4-38)  (t-g,g  +  9)r>0  VtG0. 

We  see  that  for  almost  all  values  of  A,  necessary  conditions  for  an  optimum  are 
that  (4.29)-(4.30)  and  (4.36)-(4.38)  are  satisfied.  Again,  the  system  formed  by  these 
equations  will  be  called  an  optimality  system. 

We  now  specialize  to  the  case  0  =  L2(r).  Note  that  the  hypothesis  (H10)  is 
satisfied.  Then,  using  Theorem  2.5,  we  see  that  the  inequality  (4.38)  becomes  an 
equality  and,  by  letting  z  =  t  -  g  vary  arbitrarily  in  L2(r),  we  now  have,  instead  of 
(4-38), 

(4.39)  (rg  +  |)r  =  0  VzGL2(r). 

Thus,  according  to  that  theorem,  we  have  that  for  almost  all  A,  an  optimality  system 
of  equations  is  now  given  by  (4.29)-(4.30),  (4.36)-(4.37),  and  (4.39).  However,  we  can 
go  further  and  verify  that  the  hypothesis  (Hli)  is  valid,  which  in  turn  will  justify  the 
existence  of  a  Lagrange  multiplier  satisfying  the  optimality  system  for  all  A  6  A. 

To  verify  (H 11),  we  first  note  that,  through  the  change  of  variable  f  =  T*t>,  that 
assumption  can  be  equivalently  stated  as  follows: 

ift  €  V*  satisfies  (/  +  A T*[W'(u)]*)£  =  0  and  K't  =  0,  then  £  =  0 . 

To  verify  this  version  of  (H 11),  we  assume  that  (67’=  H’(f2)  satisfies  (/  + 
A  T*  [Af'(^»)]*)4  =  0  and  K't  =  0,  i.e., 

«(Ci.Ci)  _  &(Ci,&)  +  A  ((3V>i  +  V’DCiiCi) 

+A  ((2^2)6, Ci)  =  0  VCi  €//'(«), 

«(C 2,Ca)  +  b( C21C1)  +  A  ((3^2  +  C2) 

+A  ((2^2)6. C2)  =  0  VC 2GH1(fi), 

and 

C  =  0  on  r. 

(Note  that  A'*£  =  (|r  )  Let  fi*  be  a  smooth  extension  of  ft  such  that  is  a  compact 
subset  of  i;  .  We  then  define  £ ,  if>'  and  A'  to  be  the  extension,  by  zero  outside  fi, 
of  (,  and  A,  respectively.  Let  the  forms  a'(-,  •),  b  (•,■),  and  (•,  •)'  defined  over  fi 
be  the  analogues  of  corresponding  forms  defined  over  fi.  We  may  show  from  the  last 
three  equations  that 

('gh1^'),  *GL6(n'), 

«'(Ci ,  C', )  -  6' (Ci ,  C'2)  +  A  ((3V>j2  +  V>'22)C'i  ,  Cl)' 

+A((2V>/,^,Ci)'  =  0  VCiGH’(fi'), 
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and 

-'«>.£)  +  +  A  ((3i/>iJ  +  *)€i,Ca)' 

+A  (VMNt ,  G )'  =  0  V  Ca  6  Hj(fl') . 


In  the  sense  of  distribution,  ^  satisfies 


(4.40) 


and 


-Atf  -  2A'  ■  V^j  +  (|Af  +  A  (3V>i 2  +  V-22)  -  I)*! 

-(V  •  A'  -  2A  tf'jV’ikj  =  0  in  «' 


(4.41) 


-A£  +  2A'  •  Vtf  +  (V  •  A'  +  2A  sj>\  tf&K} 

+  (|A'|2  +  A  (3^2  +  V>;2)  -  l)fa  =  0  in  0' 


We  now  quote  the  following  unique  continuation  result  whose  proof  can  be  found  in 
[17].  See  also  [12]  and  [19]. 

Lemma  4.5.  Let  Si'  be  an  open  and  connected  subset  of  Rd,  d  =  2  or  3.  Let  the 
functions  V  €  [L?oc(0')]dx,i  for  some  q  >  2  and  W  €  [Lj£~1(«')]<,xdx<l  be  given.  If 

t  €  Hloc(a') , -Afc  +  Ejul  zi=l  Wi ridb/d*i)  +  zU  Vi’t>  = 0  (» the  sense  °f 
distributions),  i  =  1, ...,  d,  and  (  =  0  on  «»  open,  non-empty  subset  of  12  ,  then  £  =  0 
on  Si' .  D 

Since  A  G  Hl(12)  and  4  G  Hl(12),  it  is  easy  to  see  that  the  coefficients  in  (4.40)- 
(4.41)  satisfy  the  regularity  requirements  of  Lemma  4.5.  Also  note  that  ^  =  0  on 
(Si  \  12)  which  contains  an  open  set.  Thus  we  obtain  that  =  0  in  12  ,  or  £  =  0  in  12, 
i.e.,  (Hll)  is  verified. 

Hence  we  conclude  that  for  all  A,  the  optimality  system  (4.29)-(4.30),  (4.36)- 
(4.37),  and  (4.39)  has  a  solution.  Thus,  we  have  Theorem  2.6  which,  in  the  present 
context,  is  given  as  follows. 

THEOREM  4.6.  Let  (^,g)  G  H*(12)  x  LJ(r)  denote  an  optimal  solution  that  minimizes 
(4.32)  subject  to  (4.29)-(4.30).  Then,  for  all  A  G  A,  there  exists  a  nonzero  Lagrange 
multiplier  if  G  H'(12)  satisfying  the  Euler  equations  (4.36)-(4.37)  and  (4.39).  □ 


4-2.3.  Verification  of  the  hypotheses  for  approximations  and  error  estimates. 

We  finally  verify  the  hypotheses  (H13)-(H19)  that  are  used  in  connection  with 
approximations  and  to  derive  error  estimates. 

A  finite  element  discretization  of  the  optimality  system  (4.29)-(4.30),  (4.36)- 
(4.37),  and  (4.39)  is  defined  in  the  usual  manner.  We  first  choose  families  of  finite 
dimensional  subspaces  Xh  C  H'(O)  and  Gh  C  L2(r)  parameterized  by  a  parameter  /» 
that  tends  to  zero  and  satisfying  the  following  approximation  properties:  there  exists 
a  constant  C  and  an  integer  r  such  that 


(4.42) 


and 


(4.43) 


inf  \\4-+h\U  <  CW|WU+i ,  V  *  G  H"*+1(12) ,  1  <  m  <  r 
4  €*» 


inf  jjz  -  zh||0,r  <  Chm 


inf 


llvllm+l/2  , 


v€H™+,/J(n),v|r=i 

V  z  €  Hm+1/2(12)|r ,  1  <m<r. 
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One  may  consult,  e.g.,  [8]  and  [15],  for  some  finite  element  spaces  satisfying  (4.42)  and 

(4.43) .  For  example,  one  may  choose  Xh  =  Vk  x  Vh  where  V*  is  the  piecewise  linear 
or  quadratic  finite  element  space  defined  with  respect  to  a  family  of  triangulations 
of  0.  In  this  case,  h  is  a  measure  of  the  grid  size.  For  simplicity  we  may  choose 
Gh  =  (Afk)|r.  i.e.,  the  functions  in  Gh  are  the  restrictions  to  the  boundary  T  of 
functions  belonging  to  Xk. 

Once  the  approximating  spaces  have  been  chosen,  we  may  formulate  the  approx¬ 
imate  problem  for  the  optimality  system  (4.29)-(4.30),  (4.36)-(4.37),  and  (4.39):  seek 
**  €  Xh,  g*  €  GK ,  and  t/h  €  Xh  such  that 

(4.44)  a(V?,^)  +  6(^,^)  +  A{((^)2  +  (^)J)^,^}  =  A(^,^)r  V  €  Vk  , 

(4.45)  a(^,^)-i(^,^)  +  A{((^)J  +  (^)2)^,^}=A(^,^)r  V*$G  Vh , 

«(tf  ,*lf)  -  +  a((3(^?)2  +  (^)V.Cf) 

+  A((2*f*$).£,tf )  =  a(*?  -  *io, C?)  V  c?  e  vk , 

«(C2h,  ri)  +  *(C2\  9?)  +  A((3(*k)2  +  (*?)»).£,<$) 

+  A((2V»?V-$hf,C2)  =  a(V>2  - *20 * C2 )  v c2h  e  vk , 


(*\g*+i|k)r  =  0  Vz*eGk. 

The  operator  Th  G  £(V ;  Xk)  is  defined  as  follows:  for  f  G  Y ,  Tkf  =  /eA*  is 
the  solution  for 

«(*!,*!  )  +  *(*a,*i)  =  </>.*?)  V*?GVk 

and 

«(*$ ,  *5)  -  *(*? ,  *$)  =  (/a ,  *2 >  V  G  Kk  . 

Since  T  =  T*.  we  define  (T*)A  =  Tk. 

We  define  the  operator  £k  :  L2(r)  — *  Gh  as  the  L2(r)-projection  on  Gh,  i.e.,  for 
each  g  G  L2(r), 

(^g,zh)r  =  (g)z',)r  Vz‘GCk. 

Since  G  =  L2(r)  is  reflexive,  Ek  is  in  fact  sin  operator  from  G*  — *  Gh. 

By  results  concerning  the  approximation  of  the  Ginzburg-Landau  equations  (see, 
e.g.,  [11]),  we  obtain 

||(T-I*)f||x-»0 

as  h  — *  0,  for  all  f  G  Y.  This  is  simply  a  restatement  of  (H13). 

(H14)  follows  trivially  from  (H13)  and  the  fsict  that  T  is  self-adjoint  and  we  have 
chosen  {T*)h  =Th. 

(H15)  follows  from  the  best  approximation  property  of  L2(r)-projections  and 
(4.43). 

(H16)  and  (H17)  follow  from  the  fact  that  N  and  T  are  polynomials.  Here  we 
also  use  imbedding  theorems  and  Cauchy  inequalities. 


(4.46) 

(4.47) 

and 

(4.48) 
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Setting  Z  =  Z  =  H1/J+‘(fi),  we  have  that  Z  (H‘(n)]’  =  AT.  For  each 
9  €  H}(0)  and  (  E  H^Q),  Sobolev  imbedding  theorems  imply  that 


[*W  n  = 


/  (3^?  +  +  (2V>1^2)»?2 

\  (3V*2  +  V»?)»ja  +  (24>ifo)m 


c  z, 


(mvr  •<)•*= 

and 


/  (fMCl  +  2t/'2C2)»?l  +  (2V»lC2)»/2  +  (2Cl  v>2)»72  \  _  T  4/3/ry,  r  7 

\(6tfaC2  +  2^1  Cl )r/2  +  (2V’lC2),7l  +  (2Cl^2)»?l  / 

(r"(4)  •  C)  •  9  =  (  *  £ )  e  L4/3(n)  c  z . 


These  relations  verify  (HI 8). 

From  the  definition  of  the  operator  K  we  see  that  K  maps  L2(r)  into  [H1^2+e(fi)]’ , 
i.e.,  K  maps  G  into  Z.  Thus  (H19)  is  verified. 

Hence,  we  are  now  in  a  position  to  apply  Theorem  3.5  to  derive  error  estimates 
for  the  approximate  solutions  of  the  optimality  system  (4.29)-(4.30),  (4.36)-(4.37),  and 
(4.39).  It  should  be  noted  that  Lemma  3.4  implies  that  for  almost  all  values  of  A,  the 
solutions  of  the  optimality  system  are  regular. 

THEOREM  4.7.  Assume  that  A  is  a  compact  interval  of  R  +  and  that  there  exists  a 
branch  {(A,^(A),g(A),9(A))  :  A  €  A}  of  regular  solutions  of  the  optimality  system 
(4.29)-(4.30),  (4.36)-(4.37),  and  (4.39).  Assume  that  the  finite  element  spaces  Xk  and 
Gk  satisfy  the  hypotheses  (4.42)-(4.43).  Then,  there  exists  a  6  >  0  and  a  h0  >  0  such 
that  for  h  <  h0,  the  discrete  optimality  system  (4.44)-(4.48)  has  a  unique  branch  of 
solutions  {(A, ^fc(A),  jjrA(A), i|fc(A))  :  A  €  A}  satisfying 


{||*h(A)  -  ^(A)|U  +  ||gk(A)  -  g(A)||0,r  +  ||9*(A)  -  9(A)||,}  <  6  for  all  A  E  A. 


Moreover, 

Iim{||**(A)  -  *(A)||i  +  ||g*(A)  -  g(A)||0>r  +  ||9*(A)  -  «KA)||i}  =  0 , 
uniformly  in  A  €  A. 

If,  in  addition,  the  solution  of  the  optimality  system  satisfies  (lft(A),g(A),9(A))  E 
Hm+1(Q)  x  Hm+1/J(ft)|r  X  Hm+1(fl)  for  A  E  A,  then  there  exists  a  constant  C, 
independent  of  h,  such  that 

IhKA)  -  *k(A)||i  +  ||g(A)  -  gfc(A)||0>r  +  119(A)  -  9h(A)||, 

<  Chm(||^(A)||m+i  +  inf  ||v||m+1/2  +  ||9(A)||m+,), 

veHm+>/»(n),v|r=g 


uniformly  in  A  E  A. 

Proof:  All  results  follow  from  Theorem  3.5.  For  the  last  result,  we  also  use  (3.25)  and 
the  estimates  (see  [11]) 


IK^T"1  -  /)*||,  <  C7»m|M|m+,  for  *  G  Hm+l(Q) , 
||((T*)k(T*)~1  -  /)9||.  =  ||(TfcT_1  -  7)911,  <  C'h"*||9||m+,  for  9  e  H”*+‘(«) 


and  (see,  e.g.,  [2],  [8]  and  [18]) 

I !(£*£-'  -  Ogllo.r  <  CbT  inf  ||v||m+1/J  for  g  €  H"*+‘/’(n)|r . 

v€H",+,/J(n),v|rag 

In  these  estimates,  the  constant  C  is  independent  of  h,  g,  t),  and  A.  □ 

Remark.  In  fact,  we  obtain  from  (4.39)  that  g  =  -q|p  which  implies 


„  inf  IMIm+l/J  <  IMIm+l/S  <  IMIm+l 

v€Hl*+,',(0),v|r=g 


so  that  the  term  (infv€H*+>/»(n),v|r=g  IMU+i/a)  in  the  right-hand  side  of  the  error 
estimate  is  redundant.  □ 

Remark.  By  using  (4.39)  again,  along  with  (4.48)  and  the  error  estimate  in  Theorem 
4.7,  we  have  the  following  improved  error  estimate  for  the  approximation  of  the  control 

g: 

llg(A)  -  gh(A)||1/2,r  <  C||g(A)  -  *ffc(A)||,  <  C/.m(||*(A)||m+1  +  |WA)||m+1) . 

Of  course,  we  also  use  the  fact  that  we  have  chosen  Gh  =  (XK)|r  C  Hl/2(r).  D 

4.3.  Dirichlet  boundary  control  for  the  Navier-Stolces  equations  of  incom¬ 
pressible,  viscous  flow 

For  this  application  we  will  use  Dirichlet  boundary  controls,  i.e.,  control  is  effected 
through  the  data  in  a  Dirichlet  boundary  condition.  Let  0  denote  a  bounded  domain 
in  R**,  d  =  2  or  3  with  a  boundary  denoted  by  T.  Let  u  and  p  denote  the  velocity  and 
pressure  fields  in  fi.  The  Navier-Stokes  equations  for  a  viscous,  incompressible  flow 
are  given  by  (see,  e.g.,  [13],  [14],  or  [20]) 

-i/ V  •  ((Vu)  +  (Vu)T)  +  (u  •  V)u  +  Vp  =  f  in  fi , 


V  ■  u  =  0  in  ft 


and 

u  =  b  +  g  on  T, 

where  f  is  a  given  body  force,  b  and  g  are  boundary  velocity  data  with  fr  b  ■  n  dT  =  0 
and  frg  ■  ndT  =  0,  and  v  denotes  the  (constant)  kinematic  viscosity.  We  have 
absorbed  the  constant  density  into  the  pressure  and  the  body  force.  If  the  variables  in 
these  equations  are  nondimensionalized,  then  v  is  simply  the  inverse  of  the  Reynolds 
number  Re. 

Setting  A  =  \/u  =  Re  and  replacing  p  with  p/A,  b  with  Ab,  and  g  with  Ag,  we 
may  write  the  Navier-Stokes  equations  in  the  form 


(4.49)  —  V  •  ((Vu)  +  (Vu)T)  +  Vp  +  Au  •  Vu  =  A f  in  fi, 


(4.50)  V  ■  u  =  0  in  fl , 
and 

(4.51)  u  =  A(b  +  g)  on  T. 
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and 


We  introduce  the  subspaces 

£§(«)  =  (peiJ(0)  |  ^pdfi  =  o| 

Hi(r)  =  {g6H'(r)  |  ^g-ndr  =  o}  . 


We  also  introduce  the  bilinear  forms 


and 


<*(u, v)  =  \  jf  ((Vu)  +  (Vu)T)  :  ((Vv)  +  (Vv)T)  d(l  V  u,  v  €  H*(0) 

6(v,«)  =  -/  qVvdd  Vv6H‘(n)andVp€lJ(n) 

Jn 


and  the  trilinear  form 


c(u,v,w)=  /(u-V)vwdft  V  u,  v,  w  e  H^ft) . 

Jn 

These  forms  are  continuous  over  the  spaces  of  definition  indicated  above.  Moreover, 
we  have  the  coercivity  properties 


(4.52) 
and 

(4.53) 


a(v, v)  >  C«,||vJ|j  Vv6H‘(D) 


sup  ^11*11°  V,€L?(n) 

0*v€HJ(O)  llvlh 


for  some  constants  Ca  and  Cj  >  0.  For  details  concerning  the  notation  employed 
and/or  for  (4.52)-(4.53),  one  may  consult  [13],  [14],  and  [20]. 

We  recast  the  Navier-Stokes  equations  (4.49)-(4.51)  into  the  following  particular 
weak  form  (see,  e.g.,  [15]):  seek  (u,p,t)  €  H^fi)  x  Ljj(fl)  x  H_1/2(r)  such  that 

(4.54)  a(u,  v)  +  i(v,  p)  -  (t,  v)r  +  A  c(u,  u,  v)  =  A  (f ,  v)  V  v  €  H1  (fi) , 


(4.55)  6(u,g)  =  0  V?€ig(n), 

and 

(4.56)  (s,u)r-A(s,g>r  =  A(s,b>r  VseH'^r). 

Formally  we  have 

t  =  [-pn  +  (Vu  +  (Vu)T)  •  n]r , 

i.e.,  t  is  the  stress  force  on  the  boundary.  The  existence  of  a  solution  (u,p,  t)  for  the 
system  (4.54)-(4.56)  was  established  in  [15]. 
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Given  a  desired  velocity  field  uo,  we  define  for  any  (u,p,t)  €  HJ(ft)  x  Lq(Q)  x 
H-‘/J(r)  and  g  G  Hj;(r)  the  functional 

(4.57)  J( u, P.  t, g)  =  ^  J  |u  -  u0|4  d(l  +  |  /r(!V.g|2  +  |g|2)  dr , 
where  V,  denotes  the  surface  gradient. 

We  define  the  spaces  X  =  H‘(fl)  x  L§(0)  x  H-*/a(r),  Y  =  (H1^)]*  x  L3( 0)  x 
H‘/J(r),  G  =  H/;(r),  and  Z  =  L3/2(fi)  x  {0}  x  H^r).  By  compact  imbedding  results, 
Z  is  compactly  imbedded  into  Y .  For  the  time  being,  we  assume  that  the  admissible 
set  0  for  the  control  g  is  a  closed,  convex  subset  of  G  —  H*  (r). 

We  then  consider  the  following  optimal  control  problem  associated  with  the 
Navier-Stokes  equations: 

(4.58)  min{J(u,p,t,g) :  (u,p,t)  G  X,g  G  6}  subject  to  (4.54)-(4.56) . 

We  define  the  continuous  linear  operator  T  G  C{Y\X)  as  follows:  for  each 
((,  r/,K)  G  Y,  T(C ,  >?,*)  =  (u,p,  t)  G  X  is  the  unique  solution  of 

a(u,  v)  +  6(v,  p)  -  (t,  v)r  =  ((,  v)  V  v  G  H1  (0) , 

*(«,  q)  =  (v,q)  V  q  G  L2(fi) 

and 

(s,u)r  =  (8,i*)r  VsGH-'^r). 

It  can  be  easily  verified  that  T  is  self-adjoint. 

We  define  the  (differentiable)  nonlinear  mapping  N  :  X  -*  Y  by 

(f-u  Vu\ 

b  ) 

or  equivalently 

(N( u,  p,  t),  (v,  q,  s))  =  -(f ,  v)  4-  c(u,  u,  v)  -  (s,  b)r  V  (v,  q,  s)  G  X 
and  define  K  :  H1'3^)  —  Y  by 


K  g  =  - 


or  equivalently 

(Kg,  (v, q, .))  =  -(s, g)r  V  g  G  H1/2(r)  ,V  (v,q,s)  G  X  . 
Clearly,  the  constraint  equations  (4.54)-(4.56)  can  be  expressed  as 
(u,  p,  t)  +  A  TN(u,  p,  t)  +  \  TK  g  =  0 , 
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i.e.,  in  the  form  (2.2).  With  the  obvious  definitions  for  ?{■)  and  £(■),  i.e., 


?(u,p,t)=  ^  J^\u-u0\Adn  V(u,p,t)€X 

and 

£(g)  =  \  ^(iv.gp  +  l*la) dT . 

the  functional  (4.57)  can  be  expressed  as 

P,  t.  g)  =  A  ^(u,  p,  t)  +  A  £ (g) , 

i.e.,  in  the  form  (2.3). 

We  are  now  in  a  position  to  verify,  for  the  minimization  problem  (4.58),  all  the 
hypotheses  of  §2  and  $3. 

4-S.l.  Verification  of  the  hypotheses  for  the  existence  of  optimal  solutions. 

We  first  verify  the  that  the  hypotheses  (H1)-(H6)  hold  in  the  current  setting. 
(HI)  is  obviously  satisfied  with  a  lower  bound  0. 

(H2)  holds  with  a  =  1  and  /?  =  2. 

(H3)  is  verified  with  the  choice  (u(°\p(°\ t(°\0)  €  X  x  0  where  (uW.p^)  is 
a  solution  to  the  Navier-Stolces  equations  with  Dirichlet  boundary  conditions,  and 
t(°)  =  [-p(°)n  +  (Vu<°)  +  (Vu<°))T)  •  n]r;  see,  e.g.,  (13]  or  [20]. 

In  order  to  verify  (H4),  we  assume  {g(n*}  C  0  C  H^(r)  is  a  sequence  satisfying 
g(")  — *  g  in  H1(r);  then  we  have  g<n)  -*•  g  in  H1/a(r)  so  that  limn_00(g(n\  v)r  = 
(g,  v)p  for  all  v  €  H^fi),  i.e.,  Kg^  —  Kg  in  Y.  Assume  that  the  sequence  {u^")}  C 
H1(ft)  satisfies  — »  u  in  H^fl);  then  — ♦  u  in  L4(fl)  by  the  compactness  of 

the  imbedding  H^Q)  L4(d).  Now, 

(N(vSn^),\)  =  c(u(n\u(n^,v)  =  c(u,u("\  v)  +  c(u^")  -  u,u^n\  v) 

— »  c(u,  u,  v)  +  0  =  (IV(u),  v)  as  n  — »  oo  . 

Hence,  (H4)  is  verified. 

The  verification  of  (H5)  follows  directly  from  the  observation  that  the  mappings 
(u,p,  t)  v-  .F(u,p,  t)  =  (1/4)  ||u  -  u||44(n)  and  g  ^  ^(g)  =  (1/2)  ||g||?  r  are  convex. 

To  verify  (H6),  we  combine  a  priori  estimates  obtained  from  the  constraint  equa¬ 
tions  and  the  functional.  Let  {u(k\p(k\  C  H*(0)  x  Lq(Q)  x  H~t/,2(r)  x 

H*  (r)  be  a  sequence  such  that 

(4.59)  j(uW  g(*))<C, 

(4.60)  a(u^fc\  v)  +  6(v,pW)  -  (t**\  v)r  +  Ac(u*‘\u(‘),v)  =  A  ( f,  v)  V  v  G  H^fi) , 

(4.61)  i(u(t),?)  =  0  Vq€Ll(Q), 
and 

(4.62)  (s,uW)r-A(s,g(k))r  =  A(s,b)r  Vser'^r). 
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First,  (4.59)  implies  that  (u**\  g^)  is  uniformly  bounded  in  L4(Q)  xH'(r).  For  each 
g(*),  we  may  choose  a  (wW,r(*))  €  H’(fi)  x  £;j(fl)  that  satisfies  the  Stokes  problem 

(4.63)  a(w<t>>v)  +  6(v,r<*>)  =  (f,v)  ¥v€H‘(!1), 


(4.64) 

6(w(t\g)  =  0  V  q  €  Lo(fi) 

and 

(4.65) 

=  A(g(*)  +  b)  on  T. 

Furthermore,  there  holds  the  estimate 


(4.66)  ||w<*>||,  <  C(j|f||0  +  ||b||1/2,r  +  |tg(*>||,,r)  • 

By  subtracting  (4.63)  from  (4.60)  with  v  =  also  using  (4.61)  and  (4.64), 

we  obtain 

a(u<*>  _  w<*),  u(*)  _  w(*>)  =  _A  c(u<k\  u(fc),  -  w<*>) 

(4-67)  =Ac(u(‘),u(‘)-w(‘)lU(‘)). 


Note  that 

|c(u^*\  u^*^)| 

=  l-  \j  «<‘>  •  ((V(u<*>  -  w<‘>))  +  (V(u<‘)  -  w  (*>))T)  •  u<‘)  dU 

<  C  ||(V(u(t>  -  w(i)))  +  (V(u(t)  -  ■w(t)))T||o  |®(*)|ll(n) 

<  ^  |(v(u(‘)  -  w(*))) + (v(u<*)  -  w(‘))f|i; + cA  |u(*)|;4(n) 
so  that,  using  (4.67),  we  have  that 

i  | (V(u<*>  -  w<*>))  +  (V(u<*>  -  w(*)))T||o2  <  C\  ||u<*)|[4(n)  . 

Then,  by  (4.66)  and  the  triangle  inequality,  we  nave  that 

||(Vu<‘>)  +  (Vu«)T ||o  <  C(||f||0  +  ||b||, .*2tr  +  ||g(fc)l|i,r  +  ||u(fc)|ll«(n)}  • 

Thus, 

||(Vu(k))  +  (VuW)T||o  +  ||uW||0lr 

<  H(Vu<‘>)  +  (Vu(*))T||o  +  ||b||o,r  +  ||g(t)||o,r 

<  C(||f||0  +  ||b||1/2,r  +  ||g(k)||..r  +  l|u(t)lllnn))  • 

Since  the  mapping  u  *-♦  j|Vu  +  (Vu)T||o  +  ||u||o,r  defines  a  norm  on  H1  (fl)  equivalent 
to  the  standard  H^Q^norm,  we  have  that 

l|u(t)lli  <  C{||f||0  +  ||b||1/2,r  +  Hg^ikr  +  ||u<*>||l«(n)}  , 
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and,  since  ||u(*)||t4(n)  and  ||*<*)Hi,  r  are  uniformly  bounded,  we  conclude  that  ||u^||i 
is  uniformly  bounded  as  well.  One  easily  concludes  from  (4.60)  that  _  i /2,r  is 

uniformly  bounded.  Thus  (H6)  is  verified. 

It  is  now  just  a  matter  of  citing  Theorem  2.1  to  conclude  the  existence  of  an 
optimal  solution  that  minimizes  (4.57)  subject  (4.54)-(4.56). 

THEOREM  4.8.  There  exists  a  (u,  p,  t,g)  €  Hx(fi)  x  Z.q(Q)  x  H-1/2(fl)  x  ©  such  that 
(4.57)  is  minimized  subject  to  (4.54)-(4.56).  □ 

4-3.2.  Verification  of  the  hypotheses  for  the  existence  of  Lagrange  multipliers. 

We  now  assume  (u ,p,  t,g)  is  an  optimal  solution  and  turn  to  the  verification  of 
hypotheses  (H7)-(H9). 

The  validity  of  (H7)  is  obvious. 

(H8)  holds  since  the  mapping  z  ►-*  £(g)  —  (1/2)  /r((V, g|2  +  |gj2)  dr  is  convex. 
(H9)  can  be  verified  as  follows.  For  any  (u,  p,  t)  6  X ,  the  operator  N'(n,p,  t)  : 
X  —*  Y  is  given  by 


(u • W +  v • Vu 

0 


for  all  (v,  q,  s)  €  HJ(f2)  x  Ll(Q)  x  H  ^2(r).  Thus  we  obtain  N'(n,p,  t)  •  (v,q,  s)  6 
L3/2(fl)  x  {0}  x  H^r)  =  Z. 

The  Lagrangian  is  given  by 


C(vL,p,t,g,v,<t>,r,k) 

=  k  J( u,  g)  -  {a(u,  v)  +  A  c(u ,  u,  v)  -f  b(t/,  p)  +  b( u,  <f>)  -  (r,  u)r 

-  (*.  v)r  -  A  (f ,  v)r  +  A  (t,  b)r  +  A  (r,  g)r} 

for  all  (u,p,  t,  g,*,  d>,r,  k)  €  X  x  G  x  X  x  R  =  Hx(«)  x  i2(fi)  x  H"1'2^)  x  H^(r)  x 
Hx(«)  x  L§(S1)  x  H-x/2(r)  x  R.  Note  that  in  this  form  of  the  Lagrangian,  the 
Lagrange  multiplier  ( v ,  <j>,  r)  £  X  —  Y*  so  that  we  have  already  introduced  the  change 
of  variables  indicated  between  (2.17)-(2.18)  and  (2.19)-(2.21). 

Having  verified  the  hypotheses  (H7)-(H9),  we  may  apply  Theorem  2.4  to  conclude 
that  there  exists  a  Lagrange  multiplier  (i/,  <j>,r)  6  X  =  H1  (fl)  x  L2(fi)  x  H-1^2(r) 
and  a  real  number  k  such  that 

(4.68)  {v,<j>,r)  +  A  T*  ([N\ u,  p,  t)]*  •  (v,  <b,  r)  -  k  J(UiP> t)(u,  p,  t ,  g))  =  0 

and 

(4.69)  £(u,p,t,z  v,<j>,T,k)  <  £(u,p,t,g,i/,^,T,fc)  Vz€0 

and  that  for  almost  all  values  of  A,  we  may  choose  4=1. 

Recall  that  T*  =T.  Also,  note  that  for  (u,p,  t)  €  X  =  H1(n)xL§(fi)xH_1^2(r), 
the  operator  [N'(u,p,  t)]*  :  X  — *  Y  is  given  by 

(— u  •  Vv  +  v  •  (Vu)T  \ 

0  1  V  (v,g,s)  6  X. 
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Thus,  (4.68),  with  Jb  =  1,  can  be  rewritten  as 

,4  a(w,  v)  +  A  c(w,  u,  v)  +  A  c(u,  w,  v)  +  6(w,  -  {t,  w)r 

=  A  ((u  -  u0)3,  w)  V  w  G  H1  (Q) , 

(4-71)  6(v,r)  =  0  VrGLg(n), 

and 

(4.72)  (y,v)r  =  0  Vy€H"1/2(r). 

In  the  right-hand  side  of  (4.70),  we  use  the  notation  (v3,w)  =  E?=j  («?,«*). 

Using  the  definition  of  the  Lagrangian  functional,  (4.69),  with  4=1,  .an  be 
rewritten  as 


|  (V,a,  V,z)r  +  ^  (z,  z)r  -  ^  (V,g,  V,g)r 

-^(g.g)r-A(r,z)r-i-A(r,g)r>0  VzG0. 

For  each  c  6  (0, 1)  and  each  z  €  0,  by  plugging  ez  +  (1  —  e)g  €  0  into  the  last 
inequality  we  obtain 

4  (v-g-  v»(z  ~  g))r  + f  (g. z  “  g)r  +  j  (V. (z  -  g),  V,(z  -  g))r 
+j (z-g,z-g)r-e(r,z-g)r  >  0  VzG0 
so  that,  after  dividing  by  e  >  0  and  then  letting  f  — *•  0+,  we  obtain 

(4.73)  (V,g, V,(z  —  g))r  +  (g, z  —  g)r  —  (r, z)r  >  0  VzG0. 

We  see  that  for  almost  all  values  of  A,  necessary  conditions  for  an  optimum  are 
that  (4.54)-(4.56),  (4.70)-(4.72)  and  (4.73)  are  satisfied.  Again,  the  system  formed  by 
these  equations  will  be  called  an  optimality  system. 

We  now  specialize  to  the  case  0  =  Hj,(r).  Note  that  the  hypothesis  (H10)  is 
satisfied.  Then  using  Theorem  2.5,  we  see  that  the  inequality  (4.73)  becomes  an 
equality  and,  by  letting  z  =  k  —  g  vary  arbitrarily  in  H^(r),  we  now  have,  instead  of 
(4-73), 

(4.74)  (V,g,  V,z)r  +  (g,  z)r  -  (r,  z)r  =  0  VzG0  =  H‘(r). 

Thus,  according  to  that  theorem,  we  have  that  for  almost  all  A,  an  optimality  system 
of  equations  is  now  given  by  (4.54)-(4.56),  (4.70)-(4.72)  and  (4.74).  However,  we  can 
go  further  and  verify  that  the  hypothesis  (Hll)  is  valid,  which  in  turn  will  justify  the 
existence  of  a  Lagrange  multiplier  satisfying  the  optimality  system  for  all  A  G  A. 

We  now  verify  (Hll)  which  we  again  note  can  be  equivalently  stated  as  follows: 

i ft  €  Y"  satisfies  (/  +  A  T*[Ar'(u)],)(  =  0  and  =  0,  then  4=0. 
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To  verify  this  hypothesis,  we  assume  that  ({,  <r,9)  £  Y*  =  H*(fi)  x  L2(fi)  x  H-1^J(r) 
satisfies  (/  -f  A7'*[iV/(u,p,  t)]*)  (£,  er,9)  =  (0,0,0)  and  K*{£,<r,9)  =  0,  i.e., 


a(w,{)  +  Ac(w,u,£)  +  Ac(u,w,£)  +  6(w,<r)  -  {9,  w)r  =  0  V  w  6  H!(fl) , 


and 


*«,r)  =  0  VreLS(n), 
(y,€>r  =  0  VyeH-1/2^), 

6  =  0  on  T. 


(Note  that  K*(f,<r,9)  =  9.)  Let  fi'  be  a  smooth  extension  of  fi  such  that  H  is  a 
compact  subset  of  fi  .  We  then  set  £ ,  <r'  and  u'  to  be  the  extension,  by  zero  outside 
fi,  of  <t  and  u,  respectively.  We  may  show  from  the  last  four  equations  that 


*'eH1(ft'),  <r'€Lg(0'), 


(4.75)  a'(w,(')  +  Ac^w.u',^')  +  Ac'(u',w,£')  +  6'(w,  o')  =  0  Vw€Hj(fi'), 
and 

(4.76)  6'«',r)  =  0  Vrelg(fi'), 

where  the  forms  a'(-,  -),  b‘  {-,  •)  and  c^-,  •,  •)  defined  over  fi  are  the  analogues  of  cor¬ 
responding  forms  defined  over  fi.  Using  a  unique  continuation  result  for  the  system 
(4.75)-(4.76)  that  was  established  in  [16]  or  [17],  we  obtain  =  0  and  o'  =  0  in  fi‘, 
or  (  =  0  and  <r  =  0  in  fi.  Thus  (H 11)  is  verified. 

Hence  we  conclude  that  for  all  A,  the  optimality  system  (4.54)-(4.56),  (4.70)- 
(4.72),  and  (4.74)  has  a  solution.  Thus,  we  have  Theorem  2.6  which,  in  the  present 
context,  is  given  as  follows. 

THEOREM  4.9.  Lei  (u,p,t,g)  £  H1(fi)xL§(fi)xH-1^2(r)xHi(r)  denote  an  optimal 
solution  that  minimizes  (4.57)  subject  to  (4.54)-(4.56).  Then,  for  all  A  6  A,  there 
exists  a  nonzero  Lagrange  multiplier  ( )  £  H*(fi)  x  Lg(fi)  x  H-1^2(r)  satisfying 
the  Euler  equations  (4.70)-(4.72)  and  (4.74).  0 

Note  that,  in  the  above  expression,  we  have  already  employed  hypothesis  (H12) 
which  in  the  current  context  is  easily  seen  to  be  satisfied  with  E  :G  —+G*  defined  by 

{Eg,  z)  =  jf  (V, g  •  v,z + g  •  z)  dr  v  z  €  K (n  =  g  . 

We  also  note  that  for  each  fixed  r,  (4.74)  with  g  £  H),(r)  is  equivalent  to 
(4-77)  (V,g,V,k)r  +  (g,k)r  +  Tjfk  ndr  =  (T,k)r  VkeH'(r) 
and 

(4.78)  J^gndT  =  0, 
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where  7  €  R  is  an  additional  unknown  constant  that  accounts  for  the  single  integral 
constraint  of  (4.78).  The  equivalence  can  be  shown  as  follows.  First,  an  application 
of  Lax-Milgram  Lemma  to  (4.74)  on  the  space  H^(T)  guarantees  the  existence  and 
uniqueness  of  a  solution  g  €  Hi(T)  to  (4.74);  this  solution  g  clearly  satisfies  (4.77)- 
(4.78)  with  7  =  /r(r  -  n  -  V,g  :  Vjii  -  g  •  n)  dr.  Conversely,  any  solution  (g,7)  of 
(4.77)-(4.78)  trivially  satisfies  (4.74).  Although  (4.74)  and  (4.77)-(4.78)  are  equivalent, 
the  latter  is  more  easily  discretized. 


4-8.8.  Verification  of  the  hypotheses  for  approximations  and  error  esittnaies. 

We  finally  verify  the  hypotheses  (H13)-(H19)  that  are  used  in  connection  with 
approximations  and  error  estimates. 

A  finite  element  discretization  of  the  optimality  system  (4.54)-(4.56),  (4.70)- 
(4.72),  and  (4.74)  is  defined  as  follows.  First,  one  chooses  families  of  finite  dimensional 
subspaces  Vh  C  H'(fi)  and  Sh  C  L2(fl).  These  families  are  parameterized  by  the  pa¬ 
rameter  h  that  tends  to  zero;  commonly,  this  parameter  is  chosen  to  be  some  measure 
of  the  grid  size  in  a  subdivision  of  (1  into  finite  elements.  We  let  Sq  =  Sh  fl  L2(A)  and 

v‘  =  vAnH‘(fi). 

One  may  choose  any  pair  of  subspaces  Vh  and  Sh  that  can  be  used  for  find¬ 
ing  finite  element  approximations  of  solutions  of  the  Navier-Stokes  equations.  Thus, 
concerning  these  subspaces,  we  make  the  following  standard  assumptions  which  are 
exactly  those  employed  in  well-known  finite  element  methods  for  the  Navier-Stokes 
equations.  First,  we  have  the  approximation  properties:  there  exist  an  integer  k  and 
a  constant  C,  independent  of  h,  v  and  q,  such  that 

(4.79)  inf  ||v  -  vh||,  <  CTHMU+i  V  v  G  Hm+1  (fl) ,  1  <  m  <  k 
and 

(4.80)  inf  H?  -  9a||o  <  Chm||g||m  V  q  G  //m(fl)  n  Lg(fl) ,  1  <  m  <  it . 


Next,  we  assume  the  inf-sup  condition,  or  Ladyzhenskaya-Babuska-Brezzi  condition: 
there  exists  a  constant  C,  independent  of  h,  such  that 


(4.81) 


.  ,  b(vh,qh) 

0  iiv*iii  Iloilo 


>c. 


This  condition  assures  the  stability  of  finite  element  discretizations  of  the  Navier- 
Stokes  equations.  For  thorough  discussions  of  the  approximation  properties  (4.79)- 
(4.80),  see,  e.g.,  [2]  or  [8],  and  for  like  discussions  of  the  stability  condition  (4.81),  see, 
e.g.,  [13]  or  [14].  The  latter  references  may  also  be  consulted  for  a  catalogue  of  finite 
element  subspaces  that  meet  the  requirements  of  (4.79)-(4.81). 

Next,  let  Ph  =  Vh|r,  i.e. ,  Ph  consists  of  the  restriction,  to  the  boundary  T,  of 
functions  belonging  to  Vfc.  For  all  choices  of  conforming  finite  element  spaces  Vh,  e.g., 
Lagrange  type  finite  element  spaces,  we  have  that  PA  C  H-,/,2(r).  For  the  subspaces 
pfc  _  vfc|r,  we  can  show  the  following  approximation  property:  there  exist  an  integer 
k  and  a  constant  C,  independent  of  h  and  s,  such  that 

IIs  -  •*ll-,/2.r 

(4  821  *  ^ 

1  ;  <Chm  inf  ||v||m  V  s  G  Hm(fl)|r ,  1  <  m  <  t . 

veH“(n),v|r=« 
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We  also  use  the  following  inverse  assumption:  there  exists  a  constant  C,  independent 
of  h  and  ah ,  such  that 


(4.83)  )\sk\\.x  <  C/i'-«||s%,r  V  s'1  G  P\  -1/2  <  q  <  s  <  1/2  . 

See  [2]  or  [8]  for  details  concerning  (4.82)  and  (4.83).  See  also  [15]  for  (4.82). 

Now,  let  Q*  =  Vh|r  ,  i.e.,  Qh  consists  of  the  restriction,  to  the  boundary  T,  of 
functions  belonging  to  VA.  Again,  for  all  choices  of  conforming  finite  element  spaces 
Vh  we  then  have  that  C  H^T).  We  can  show  the  approximation  property:  there 
exist  an  integer  k  and  a  constant  C ,  independent  of  h  and  k,  such  that  for  1  <  m  <  k, 
0<s<l  andkeHm+1(fi)|r, 


(4.84) 


kk€QO 


This  property  follows  from  (4.79),  once  one  notes  that  the  same  type  of  polynomials 
are  used  in  Qh  as  are  used  in  \K.  We  set  Gh  =  Qh  fl  H],(r). 

Once  the  approximating  subspaces  have  been  chosen  we  seek  uh  E  Vh ,  ph  G  Sq, 
th  G  P'*,  gh  e  Q\  vh  G  Vh,  4>h  G  Tk  e  P\  and  7*  G  R  such  that 

(4.85)  a(uh,  vh)  +  Ac(ufc,uh,  vk)  +  b(yh,ph)  —  (vfc,  th)r  =  A  (f ,  vA)  V  vh  6  Vh  , 


(4.86) 

6(u\,h)=0  Vqh6S5, 

(4.87) 

(uh,sh)r  -  \(gh,ah)r  =  A(b,  sh)r  V  s'1  G  P\ 

(4.88) 

(V,g\V,kh)r  +  (g\kh)r  +  7'‘  ^k'*ndr=(Th,kh)r  Vkh  gQ\ 

(4.89) 

J^ghndT  =  Q} 

(4.90) 

a('wh,vh)  +  Xc(wh,uh,vh)  4-  \c(uh,wh,vh)  +  b(wh,<f>h)  -  (w',,rA)r 

=  A  (( uh  -  uo)3,  wh)  V  wh  6  Vh  , 

(4.91) 

6(i/\rh)  =  0  V  rh  €  Sq  , 

and 

(4.92) 

(vk,yh)=0  Vy*eP*. 

Note  that  if  (4.85)-(4.92)  are  satisfied,  then  necessarily  gh  E  Gh.  Also,  in  the  right- 
hand  side  of  (4.90),  we  use  a  notation  similar  to  that  used  in  the  right-hand  side  of 
(4.70). 
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The  operator  Tk  G  £(V',  Xk)  is  defined  as  the  solution  operator  for 
a(u* ,  v*)  +  b(vk ,  pk)  -  <vfc ,  tk)r  =  <f ,  vk)  Vv‘eVk, 

6(u\,ft)  =  0  VqheSk, 
and 

{u*,s*)r  =  {b,s*)r  VsfcGP\ 

i.e,  for  each  f  €  Y,  Tht  =  **  e  Xh  is  the  solution  of  the  above  system  of  equations. 
Since  T  -  T,  we  define  (T*)fc  =  Tk. 

We  define  the  operator  Ek  :  G*  — »  Gh  as  follows.  For  each  rgG’,  g*  =  £Ar  if 
and  only  if 


(V,g\V,*',)r  +  (g'i,*k)r  +  7 k  JT*h  nc/r  =  (r\*A}r 


V**€Q'1 


and 


•  n  dr  =  0 . 


The  existence  and  uniqueness  of  a  solution  (gk,yk)  6  Qk  x  R  is  guaranteed  by  the 
Brezzi  theory  for  mixed  finite  element  methods  (see  [4]  or  [5])  and  the  inequalities 


(4.93)  (V,k*,  V,k*)r  +  (kfc, kh)r  >  C|jkA|Jjr  V  kfc  G  Q*  C  H1  (r) 

and 


(4.94) 


sup 

o^k»gq» 


yk  frkk  -n  dr 
||k*||i,r 


>  cirhl 


Vr*€R. 


The  solution  necessarily  satisfies  gk  G  Gk.  Thus  the  operator  Eh  is  well  defined. 

With  these  definitions  we  see  that  (4.85)-(4.92)  can  be  written  in  the  form  (3.1)- 
(3.3). 

By  results  concerning  the  approximation  of  the  Navier-Stokes  equations  with 
inhomogeneous  boundary  conditions  (see  [15]),  we  obtain 

||(T  —  T*)/||x  — » 0 


as  h  — ►  0,  for  all  /  =  ((,  i),  it)  G  Y.  This  is  simply  a  restatement  of  (H13). 

(H14)  follows  trivially  from  (H13),  the  fact  that  T  is  self-adjoint,  and  the  choice 

(y*)h  _  rph 

To  verify  (H15),  we  note  that  the  nondiscretized  version  of  (4.93)-(4.94)  certainly 
also  holds,  i.e., 


and 


(V.k,  V,k)r  +  (k,  k)r  >  C||k||i  r  V  k  G  H1  (r) 


8Up  ~^l  V 1  6  R  • 

o*k6H>(r)  ||k||i,r 


Using  the  Brezzi  theory  for  mixed  finite  element  method  (see  [4]  or  [5]),  we  obtain 
that 

||(JS7  —  £*)r||llr -♦  0  asA-»0, 
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which  verifies  (H15). 

(H16)  and  (H17)  follow  from  the  fact  that  N  and  T  are  polynomials.  Here  we 
also  use  imbedding  theorems  and  Cauchy  inequalities. 

We  set  Z  —  L3/2(fl)  x  {0}  x  {0}.  For  each  (v,$,s)  6  X  =  H*(Q)  x  L§(f2)  x 
H~1/2(r)  and  (w,r,k)  6  X  =  H1^)  x  L2(Q)  x  H"1/2(r),  Sobolev  imbedding  theo- 
rems  imply  that 

(-(uV)v  +  v(Vu)T\ 

0  I  €  Z, 

/  -(w  •  V)v  +  v  •  (Vw)T  \ 

([N"(u,p,t)]m  (v,q,s))  •  (w,r,k)  =  -  I  0  j€Z, 

and 

3(tll  -Uoi)2^!!)! 

3(u<j  -  uojywr'd 
0 

0 

where  d  (=  2  or  3)  is  the  space  dimension.  These  relations  verify  (H18). 

From  the  definition  of  the  operator  if  we  see  that  K  maps  H*  (r)  into  L3/,2(Q)  x 
{0}  x  H1^),  i.e.,  K  maps  G  into  Z.  Thus  (H19)  is  verified. 

Hence,  we  are  now  in  a  position  to  apply  Theorem  3.5  to  derive  error  estimates 
for  the  approximate  solutions  of  the  optimality  system  (4.54)-(4.56),  (4.70)-(4.72)  and 
(4.74).  It  should  be  noted  that  Lemma  3.4  implies  that  for  almost  all  values  of  A,  the 
solutions  of  the  optimality  system  are  regular. 

THEOREM  4.10.  Assume  that  A  is  a  compact  interval  of  R  +  and  that  there  exists 
a  branch  {(A,  u(A),  p(A),  t(A),g(A)KA),  $(A),  r(A))  :  A  £  A }  of  regular  solutions  of 
the  optimality  system  (4.54)-(4.56),  (4.70)-(4.72),  and  (4.74).  Assume  that  the  finite 
element  spaces  and  Gk  satisfy  the  hypotheses  (4.79)-(4.84).  Then,  there  exists  a 
6  >  0  and  an  h0  >  0  such  that  for  h  <  ho,  the  discrete  optimality  system  (4.85)-(4.92) 
has  a  unique  branch  of  solutions  {(X,nh{\),ph(X),th(X),gh(X),l/h(X),4>h{X),Th(X))  : 
A  €  A)  satisfying 

(||u(A)  -  uh(A)||i  +  ||p(A)  -  ph(A)||0  +  ||t(A)  -  t*(A)||_1/3ir 
+  ||g(A)  -  gA(A)||i,r  +  IKA)  -  *a(A)||,  +  ||*(A)  -  ^(A)||0 
+  ||r(A)  -  rA(A)||_I>2,r)  <  S  for  all  A  €  A. 

Moreover, 

lim  (||u(A)  -  ufc(A)||,  +  ||p(A)  -  pk(A)||0  +  ||t(A)  -  tfc(A)||_1/2,r  +  ||g(A)  -  gh(A)||1,r 
+  |KA)  -  i^(A)||,  +  ||*(A)  -  ^a(A)||o  +  MX)  -  rh(A)||_I,2,r)  =  0 
uniformly  in  A  €  A. 
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//,  in  addition,  the  solution  satisfies  (u(A),p(A),  t(A),g(A),y(A),  ^(A),r(A))  £ 
Hm+1(fl)  x  Hm{to)  x  Hm(fl)|r  x  Hm+1(n)|r  x  x  Hm(Q)  x  Hm(Q)|r  for 

A  £  A,  then  there  exists  a  constant  C,  independent  of  h,  such  that 


( ||u(A)  -  u*(A)||,  +  \\p{X)  -  p*(A)||o  +  ||t(A)  -  tfc(A)||_1/2,r  +  ||g(A)  -  gk(\)\U,r 

+  HA)  -  **(A)||,  +  ||0(A)  -  0*(A)||o  +  MX)  -  r*(A)||_1>2,r) 

<  CAm-1^2(||u(A)||m+i  +  ||p(A)||m  +  cn  inf  ||v||m 

'  v6Hm(n),v|r=t 


+  C„m+}?JL  llvl|m+l  +  |HA)l|n,+l  +  ||^(A)||m  + 
vgH™+l(n),v|r=g 


inf 

w€H»(n),w|r=r 


uniformly  in  A  £  A. 

Proof:  All  results  follow  from  Theorem  3.5.  For  the  last  result,  we  also  use  (3.25)  and 
the  estimates  (see,  e.g.,  [16]  or  [17]) 

\\(ThT~l  -  /)(u,p,t)||x  <  Chm(IM|m+1  +  ||p||m  +  £Hminf  ||v||m) 

v€Hm(0),v|r=t 

for  u  £  Hm+1(n),  p  £  Hm(Sl),  and  t  6  Hm(n)|r , 


||((T*)'‘(T*)-1  -  /)(*,*, r)||y.  =  ||(TkT-‘  -  I)(*,4>,t)\\x 

<  CAm(||vj|fn+i  +  ||^||m  +  c  inf  |H|m) 

w6Hm(n),w|r=f 

for  v  £  Hm+1(n),  <t>  G  Hm((l),  and  r  £  Hm(ft)|r  , 

and 

\\(EhE-'-I)g\Ux<Ch"'-1'*  inf  ||v||m+i  for  g  £  Hm+1  (n)|r  ■ 

v6Hm+,(0),v|r=g 

In  these  estimates,  the  constant  C  is  independent  of  h,  u,  p,  t,  g,  v,  <f>,  r,  and  A.  □ 
Remark.  If  the  control  g  £  Hm+3/2(Q)|r,  then  the  exponent  of  h  in  the  error  estimate 
of  Theorem  4.10  can  be  increased  from  (m  —  1/2)  to  m.  □ 


5.  Conclusions 

We  have  set  up  an  abstract  framework  for  the  analysis  and  approximation  of 
a  class  of  nonlinear  optimal  control  and  optimization  problems.  Nonlinearities  can 
occur  in  both  the  objective  functional  and  in  the  constraints.  Within  the  framework  we 
have  defined  an  abstract  nonlinear  optimization  problem  posed  on  infinite  dimensional 
spaces,  defined  an  approximate  problem  posed  on  finite  dimensional  spaces,  and  listed 
a  number  of  hypotheses  concerning  the  two  problems.  We  then  have  shown  that 
optimal  solutions  exist  and  that  Lagrange  multipliers  may  be  used  to  enforce  the 
constraints.  We  then  used  the  Lagrange  multiplier  rule  to  derive  an  optimality  system 
from  which  optimal  states  and  controls  may  be  deduced.  We  then  derived  existence 
results  and  error  estimates  for  solutions  of  the  approximate  problem.  The  abstract 
framework  and  the  results  derived  from  that  framework  were  then  applied  to  three 
concrete  control  or  optimization  problems  and  their  approximation  by  finite  element 
methods.  The  first  involves  the  von  Karman  plate  equations  of  nonlinear  elasticity, 
the  second  the  Ginzburg-Landau  equations  of  superconductivity,  and  the  third  the 
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Navier-Stokes  equations  for  incompressible,  viscous  flows.  It  is  certainly  possible  to 
apply  the  abstract  results  that  we  have  derived  to  a  variety  of  optimal  control  problems 
arising  in  other  settings. 
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